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1  3.  ABSTRACT 

The  objective  of  the  program  was  to  evaluate  feasibility  of  non-destructive 
electronic  measurement  (ard  logging]  of  rock  stress  in  situ,  utilizing  atomic 
resonances  that  are  pressure-sensitive.  Nuclear  quadrupole  resonance  (NQR) 
appeared  clearly  applicable,  since  no  magnetic  aligning  field  is  needed,  there 
are  several  common  NQR-active  nuclear  species  in  major  rock-forming  minerals, 
and  rf  magnetic  fields  between  0.3  and  3  megahertz  (the  resonant  frequencies) 
propagate  well  through  many  rocks.  There  are  also  alternative  resonance  methods 
which  may  be  applicable,  including  NMR,  EPR,  and  EFR.  A  comprehensive  literature 
search  has  been  completed  on  NQR  in  minerals,  which  shows  that  aluminum-27  and 
a  few  other  elements  in  feldspars  are  promising  for  NQR  measurement  of  rock 
stress.  Efforts  to  detect  these  resonances  in  mineral  specimens  were  mostly  un¬ 
successful  because  of  their  extreme  weakness.  The  outlook  for  NQR  is  now  con¬ 
sidered  marginal.  Short  studies  of  other  rf  spectrographic  methods  showed  that 
ferromagnetic  effects  cannot  be  used,  but  that  electron  paramagnetic  resonance 
(EPR)  holds  some  promise  as  a  means  of  measuring  rock  stress,  since  strong  EPR 
signals  are  readily  obtained  from  most  rocks. 
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SUMMARY  ABSTRACT 


A  desk  and  laboratory  study  has  been  undertaken,  to  determine  the  feasibility 
of  measuring  (and  logging)  rock  stress  in  situ,  by  nondestructive  electronic 
means  employing  atomic  resonances  sensitive  to  external  pressure.  A  compre¬ 
hensive  literature  search  and  study  showed  that  a  large  amount  of  research 
had  been  done  on  electronic  (radio  frequency)  spectroscopic  methods,  particu¬ 
larly  NMR,  EPR  and  "pure"  NQR.  Nevertheless,  most  of  this  work  had  been  de¬ 
voted  to  synthetic  and  organic  substances;  many  of  the  facts  required  on 
mineral  resonances  for  rock  stress  measurement  had  not  been  gathered.  An  ex¬ 
tensive  experimental  effort  was  performed,  both  in  company  laboratories  and 
in  outside  institutions.  It  has  now  been  shown  that  the  only  NQR  signals  of 
possible  value  are  those  of  aluminum,  sodium,  and  potassium  in  the  feldspars; 
but  these  are  so  weak  that  their  applicability  to  rock  stress  measurement  is 
in  doubt.  Resonances  of  ferromagnetic  minerals  are  not  usable.  The  NMR 
process  still  has  possibilities,  but  will  be  difficult  to  apply.  However, 
EPR  effects  are  stron6  and  widespread  in  minerals.  The  outlook  for  their 
eventual  application  to  measuring  rock  stress  is  considered  reasonably  good. 
Further  investigation  of  the  EPR  method  is  recommended. 
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ELECTRONIC  MEASUREMENT  OF  ROCK  STRESS 


1.  INTRODUCTION 


The  purpose  of  this  study  was  to  determine  the  feasibility  of  measuring  rock 
stress  in  situ  by  purely  electronic,  non-destructive  means.  Such  a  cay  ability 
might  permit  rapid  measuring  of  the  state  of  stress  in  the  rock  at  the  exposed 
surfaces  of  an  underground  shaft,  drift  or  tunnel,  as  well  as  continuous  log¬ 
ging  of  rock  stress  in  deep  drill  holes. 

The  phenomena  to  be  exploited  were  atomic  magnetic  resonances  at  radio  fre¬ 
quencies.  RE  fields  introduced  at  a  rock  surface  can  penetrate  to  substantial 
depths,  of  the  order  of  meters  to  hundreds  of  meters  in  most  rocks  (see 
Appendix  1).  They  may  thereby  couple  electromagnetically  with  multitudes  of 
resonant  particles  in  the  rock,  which  comprise  natural  stress  transducers. 
The  variety  of  atomic  resonance  of  greatest  potential  value  for  this  purpose 
was  considered  to  be  "pure  nuclear  quadrupole  resonance"  or  NQR,  which  is 
applicable  to  many  naturally-occurring  nuclear  species  and  does  not  require 
any  magnetic  aligning  field.  Other  possible  varieties  were  nuclear  magnetic 
resonance  (NMR) ,  electron  paramagnetic  resonance  (EPR)  and  zero-field  nuclear 
magnetic  resonance  (ZNR) ,  the  first  2  employing  the  earth's  magnetic  field 
(which  generally  has  the  required  high  degree  of  uniformity),  or  an  artificial 
field  of  up  to  a  few  hundred  oersteds,  applied  by  loops  or  solenoids  carrying 
direct  current,  to  align  the  particles.  These  types  of  resonance  are  all 
applicable  to  solid  compounds  and  are  affected  by  distortions  of  the  crystal 
lattice,  hence  by  the  absolute  mechanical  stresses.  The  resonant  frequencies, 
and/or  absorption  line  widths,  and/or  line  strengths  were  expected  to  serve 
as  indicators  of  stress  in  a  rock.  The  resonance  response  of  each  rock  type 
to  stress  would  be  determined  by  a  "calibration:"  applying  known  stresses  to 
a  core  sample  in  the  laboratory  while  measuring  the  resonances,  to  permit 
interpretation  of  the  field  measurements. 

The  history  of  the  idea  of  measuring  rock  stress  with  atomic  resonances 
apparently  began  with  Mr.  0.  K.  Kowallis,  a  former  Teledyne  employee,  who 
outlined  it  in  an  internal  company  report  (Kowallis,  1968).  One  earlier  ele¬ 
ment  of  the  idea  was  Benedek's  (1965)  article  on  pressure  effects  in  iron. 

In  the  absence  of  the  "rock  calibration"  idea  mentioned  above,  Kowallis 
reasoned  that  only  a  nuclear  effect  such  as  NMR  could  be  used.  Furthermore, 
he  needlessly  ruled  out  the  possibility  of  using  the  earth's  magnetic  field 
or  an  applied  field  for  particle  alignment,  hence  confined  his  attention  to 
NMR  of  the  ferromagnetic  materials  having  strong  internal  fields,  particularly 
compounds  of  iron,  since  ^Fe  has  a  magnetic  moment  and  has  a  Larmor  resonance 
at  what  he  considered  an  appropriate  frequency  (near  45.5  MIlz)  in  pure  iron 
granules.  Kowallis  recommended  that  the  62.1  MHz  resonance  he  calculated  for 
magnetite  be  confirmed  (recognizing  that  this  mineral  is  widely  distributed 
in  nature),  and  predicted  that  under  uniaxial  stress,  the  ^Fe  resonance  line 
would  be  split  into  two.  From  data  by  Benedek  (1965)  he  estimated  a  pressure 
efrect  in  iron  of  Af/Ap  =  D.52  Hz/psi  ,  which  should  permit  measurement  of 
pressures  by  means  of  frequency  to  a  precision  of  the  order  of  200  psi  in 
favorable  instances. 
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With  modest  company  support,  Kowallis  was  able  to  detect  the  ^7Fe  resonance 
of  pure  (carbonyl)  iron  powder  at  45.4686  MHz  with  a  volume  concentration  of 
iron  as  low  as  1%,  using  a  superregenerative  spectrometer.  He  also  detected 
a  resonance  in  an  Ottowa  sand  sample  at  45.250228  MHz,  but  could  not  find  the 
predicted  62.1  MHz  resonance  of  magnetite,  nor  any  resonances  in  several  other 
soils  and  rocks. 

Mr.  Kowallis'  report  was  subsequently  studied  at  the  Canadian  Department  of 
Energy,  Mines  and  Resources  in  Ottowa,  and  a  detailed  report  was  issued  (Gray, 
1968).  This  report  pointed  out  that  the  Ottowa  sand  'esonance  should  not  arise 
from  the  0.02%  of  Fe203  in  the  St.  Peter  sandstone  of  Illinois  and  its  deriva¬ 
tives,  including  Ottowa  sand.  Instead,  it  was  suggested  that  this  resonance 
coull  be  a  zero-field  ("pure")  nuclear  quadrupole  resonance  (NQR)  of  the 
oxygen-17  comprising  0.019%  of  the  sand,  or  that  of  the  aluminum-27  comprising 
0.03%  of  the  sand.  This  suggestion  was  based  on  the  much  stronger  magnetic 
moments  of  these  nuclei  as  well  as  their  greater  abundance  compared  to  the 
0.0003%  of  57Fe  in  Ottowa  sand;  no  comparable  data  were  available  on  the 
strength  of  NQR  signals  from  quartz  (^70)  and  corundum  (27a1)  at  the  time. 

Gray  (1968)  also  noted  that  uniaxial  stress  would  broaden  and  weaken  resonance 
lines,  rather  than  split  them,  that  each  rock  might  require  "calibration"  with 
pressure  in  the  laboratory,  and  that  borehole  signals  should  become  strongest 
when  the  hole  was  pressurized  to  the  undisturbed  stress  level  of  the  rock. 

On  the  basis  of  Kowallis'  work  and  Gray's  recommendations,  Teledyne  Geotech 
proposed  to  ARPA  in  1971  to  extend  their  research,  as  the  basis  of  a  possible 
tool  to  aid  in  the  Rapid  Excavation  program.  Because  of  a  recognized  need 
for  better  rock-stress  measuring  methods,  the  present  program  of  research  was 
supported  despite  the  expectation  by  some  in  the  Government  that  the  probability 
of  success  was  less  than  50%. 

The  program  has  been  concentrated  almost  entirely  upon  pure  nuclear  quadrupole 
resonance  (NQR) . 

Practically  all  our  results  have  been  negative;  we  conclude  that  stress  meas¬ 
urements  in  common  types  of  rock  in  the  field  by  NQR  are  not  feasible  at  this 
time.  However,  the  possibility  remains  open  with  regard  to  certain  rocks,  if 
electron  paramagnetic  resonance  (F.PR)  or  NMR  of  selected  minerals  is  employed. 
Further  research  on  these  methods  should  be  performed. 

The  following  sections  summarize  the  results  of  a  very  thorough  literature 
search  (covering  over  one  mi1 lion  papers  in  physics  and  chemistry  published 
since  1950),  and  study  (of  about  200  selected  papers),  NQR  experiments  in 
the  most  advanced  equipment  available  today,  and  other  experiments  of  a 
more  preliminary  character  i  NMR  and  RPR. 


2.  PURE  NUCL  QUADRUPOLE  RESONANCE  (NQR) 


2.1  BASIC  PRINCIPLES 

Good  summaries  of  NQR  principles  are  given  in  numerous  review  papers.  Possibly 
the  best  are  the  monograph  by  Das  and  Hahn  (1958)  and  the  recent  summary 
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chapter  by  Schempp  and  Bray  (1^701.  For  our  purposes,  a  very  simple 
qualitative  outline  will  suffice. 

The  nuclei  of  certain  atoms  behave  as  if  their  electric  charge  were  unevenly 
distributed,  over  either  a  prolate  or  oblate  spheroid  instead  of  a  sphere. 

The  symmetry  axis  of  the  spheroid  is  also  the  spin  axis,  as  shown  in  figure  1. 
The  "protruding"  parts  are  +  with  respect  to  the  intermediate  areas  of  the 

nucleus,  so  the  result  is  a  square  electric  quadrupole:  -  +  ,  which  will  ex¬ 
perience  a  torque  in  the  presence  of  an  electric-field  gradient.  Such  grad¬ 
ients  are  produced  around  the  nucleus  by  the  orbital  electrons  of  the  atom, 
and  by  charged  particles  in  other  surrounding  atoms  of  the  crysta1  lattice. 
Hence,  in  crystalline  solids  there  4 s  a  definite,  builc-in  aligning  field, 
and  the  spinning  quadrupolar  nucleus  will  precess  around  it  at  a  character¬ 
istic  Larmor  frequency.  If  an  alternating  magnetic  field  of  the  same  fre¬ 
quency  is  applied  to  the  crystal  by  an  external  coil,  enormous  numbers  of 
identical  "NQR-tuned"  nuclei  can  absorb  energy  from  it,  thus  affecting  the 
impedance  of  the  coil  in  a  measurable  way. 

The  rules  of  quantum  mechanics  govern  the  entire  process.  Only  nuclei  with 
"spin  numbers"  I  of  1  or  greater  have  an  electric  quadrupole  moment.  1  nergy 
absorption  occurs  through  quantum  transitions.  The  resonance  frequencies  v 
are  determined  in  simple  cases  by  the  magnitude  of  the  nuclear  quadrupole 
moment  eQ,  measured  in  (esu)c.m2,  the  field  gradient  eq  in  (esu)cm  ,  and  the 
quantum  numbers  I  and  m,  as  follows: 


v  =  (E-j^  -  E2)/h 


where 


Bi 


3mj2  -  1(1+1) 
4  I  (21-1) 


(1) 


(2) 


Between  the  two  energy  levels  Ej  and  E2,  the  quantum  number  m^  must  change 
by  unity.  The  expression  in  brackets  in  (1)  is  a  coefficient  that  is  char¬ 
acteristic  of  the  substance,  called  the  Nuclear  Quadrupole  Coupling  Constant 
(NQCC) .  It  is  determined  by  experiment,  or  in  favorable  cases  it  can  be 
found  by  elaborate  calculations  if  Q  is  known.  In  many  substances,  the  ten¬ 
sor  q  is  not  isotropic,  and  equation  (1)  becomes  inaccurate.  Table  1  lists 
approximate  equations  that  may  be  used  to  find  the  resonance  frequencies  v, 
providing  the  gradient-assymmetry  parameter  n  is  not  too  large  (from  0  to  0.3, 
perhaps).  Spin  numbers  I  to  be  used  for  various  nuclear  species  are  given 
in  table  2.  Values  of  the  Nuclear  Quadrupole  Coupling  Constant,  in  megahertz, 
and  of  the  parameter  i  for  some  substances  of  interest  will  be  given  later  in 
table  4. 

Since  the  field  gradient  q  depends  on  the  distances  of  the  charges  in  adjacent 
atoms,  it  is  sensitive  to  changes  in  the  volume  of  the  crystal.  Hence,  the 
NQR  frequencies  are  temperature  and  pressure-sensitive.  The  latter  is  the 
phenomenon  we  hope  to  exploit  for  the  electronic  measurement  of  absolute 
stresses  on  rocks  in  situ. 
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Figure  1.  Apparent  shapes  of  quadrupolar  nuclei 


Figure  2.  NQR  spectrometer  circuit  used  (from  Gill,  et  al ,  1967) 
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iable  1.  Transition  frequencies  between  quadrupole  energy  levels 
_ (from  Biryukov,  et  al ,  1969) _ 


Nuclear  Transi- 
spin  /  lion 


Frequencies 


1 

+1^0 

0  pi  —  1 

3 

+  1  pin 

V!-  (1  +  2*i  +  1.4583ti2-(i-0521t)»-0.7U20n4) 

-1  pin 

v2  =  -^TPT-  (*  -  2rl  +  14S83*l2-f  0.0521*1* -0.7020*1*) 
2()ft 

+2  pi  +1 

v8  =  -^f  (3-2*1  +  0.2081  |2  +  0.0521*1* +  0.O076t|*) 

1 

to 

Tl 

l 

v4  =  (3  +  2*1  +  1 .8750*)2 -  0.0521  *)»  -  0-6868*i« ) 

20ft 

+3  pi +2 

(5  +  0.2083*i2  +  0.0521*i3  +  0.0076ri«) 

20ft 

-3  pi  -2 

v«-  ~npr-  (5— 1.4583*!2— 0.0521*13  +  0.7020*14) 

2l)ft 

3/2 

3/2  pi  1/2 

5/2 

1/2  pi  3/2 

v,=  (1  +  1.0926*)2— 0.6340*1«) 

2Uft 

3/2  pi  5/2 

Vj=  (1  0.2037*12  +  0.1622*14) 

20ft 

7/2 

1/2  pi  3/2 

v,=  (1  +  3-6333*l2—  7 .2607 *|4) 

14ft 

3/2  pi  5/2 

v2=  (1  05667n2  +  1.8595*i4) 

14ft 

5/2  Pi  7/2 

v3  =  (1  U.1000*!2— 0.0180*1<) 

9/2 

1/2  pi  3/2 

V,  =  (1  +  9 .0333 *1*  —  4  5  -691 0  *l4 ) 

3/2  pi  5/2 

Vj  =  (1  -1 .3381 *l2  4- 1 1  7224*l4) 

24ft 

5/2  pi  7/2 

v3=-3:;^  (1  0.1857*1* -0.1233*1*) 

7/2  pi  9/2 

v4  =  (1  -  n  0809*1*  -  0  0043*|4) 
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Element 


1 


Table  2.  Characteristics  of  the  stable  isotopes  of  nuclei  of  spin  /  > 
_ (from  Biryukov,  et  al,  1969) 


l— 

o 

X 

E 

5 

c 

ix 

o 

X 

E 

1 

r 

X 

v"- 

c 

a, 

to 

Magnetic 

moment, 

O 

>.  • 

■~r  c/i 

IN 

—  E 

c  .  o 

o 

c 

rj  * 

^  ru 

0) 

nuclear 

2  3 
c  -  So 

3  C  3 
*-  w  1 

E 

0 

< 

«/) 

2 

3  O 

5 

x  "a 

u 

3 

z 

magnetons 

5c  O  . 

S2  “  M 

^  X 

*3  C  o 

ra  3  rM 

3  2  ' 

O  E  ^ 

Ouadrupole 

moment 

ratio 


D 

1 

■) 

00156 

1 

(1.85738 

r:»  3.5 

0.00277 

Li 

3 

f. 

7.43 

1 

0.82191 

626.5 

0.00046 

Qo/Qi  = 

7 

92  57 

3/2 

3.2559 

1654  7 

-0.042 

=  0019 

Be 

4 

9 

100 

3/2 

-1.1774 

598.3 

0-02 

B 

5 

10 

18  83 

3 

1.8006 

457  5 

0.074 

Q\ol  Q\  i  ■= 

11 

81.17 

3/2 

2.6880 

1366 

0.0355 

=  2.084 

N 

7 

14 

99.635 

1 

0  40357 

307.6 

0.02 

0 

8 

17 

0037 

5/2 

—  1 .8930 

577.2 

-0-026 

Na 

11 

23 

100 

3/2 

2.2161 

1126.2 

0.1 

Mg 

12 

25 

1005 

5/2 

— 0  85471 

2600 

Al 

13 

27 

100 

5/2 

3.6385 

1109-4 

0149 

S 

16 

33 

0.74 

3/2 

0.64274 

326.6 

-0.064 

Cl 

17 

35 

75.4 

3/2 

082089 

417.2 

—00797 

CW  < ?37  = 

37 

24.6 

3/2 

0.68329 

347  2 

— 0-0621 

—  1-2688 

K 

19 

39 

93.08 

3/2 

0-39094 

198.7 

0.09 

— 

41 

6.91 

3/2 

0.21453 

1092 

0.11 

=  1220 

Ca 

20 

43 

0.13 

7/2 

-1.3153 

286  5 

Sc 

21 

45 

100 

7/2 

4-7491 

1034.3 

Ti 

22 

47 

775 

5/2 

-0.7871 

240.0 

49 

551 

7/2 

-1.1023 

240.1 

V 

23 

51 

100 

7/2 

5.1392 

11183 

0-3 

Cr 

24 

53 

9.5 

3/2 

—0.4735 

2406 

Mn 

25 

55 

100 

5/2 

3.4610 

1055-3 

05 

Co 

27 

59 

100 

7/2 

4.6388 

1010.3 

05 

<u 

29 

63 

69  1 

3/2 

2.2206 

1128-5 

■  015 

Qnl  Qtf>  — 

65 

30  9 

3/2 

2.3790 

1209 

-0.14 

=  10806 

Ga 

31 

69 

602 

3/2 

20108 

1021.8 

,2  19 

QttIQn  — 

71 

39.8 

3/2 

2.5549 

1298  4 

0.12 

=  1.5867 

Ge 

32 

73 

761 

9/2 

-0.8768 

1485 

-0.2 

As 

33 

75 

100 

3/2 

1.4349 

729.2 

03 

Br 

35 

79 

506 

3/2 

20990 

1066.7 

033 

Qnl  Q» i  = 

81 

49.4 

3/2 

2  2626 

1149.8 

0.28 

=  1  1970 

Kr 

36 

83 

11.55 

9/2 

-0-968 

164 

0.22 

Rb 

37 

85 

72.8 

5/2 

1.3483 

411.1 

0.3 

Qn/Qtn  — 

87 

27.2 

3/2 

27415 

1393  2 

0.15 

—  2  07 

I 
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Table  2,  (Continued) 


<u 

x: 

, 

Magnetic  1 
morrent. 

u 

E 

c 

£ 

x> 

E 

3 

C 

EL 

to 

k. 

>» 

w: 

a 

•H 

1 

i/j 

V) 

o 

rg 

c 

»  ■«» 

Ouadrupole 

c 

0> 

E 

o 

E 

a 

c 

VI 

n 

3 

<U 

£ 

nuclear 

0> 

c 

be 

X 

o 

3 

nj 

bo 

a. 

X) 

?  3 

1 

p  ° 

moment 

r.mo 

w 

UJ 

0 

<  * 

nj 

2 

C 

“3  | 

2 

magnetons 

CO 

N 

nj 

3 

o 

p 

0  . 

E  O 

Sr 

38 

87 

7.02 

9/2 

-1.0893 

184.5 

Nb 

41 

93 

100 

9/2 

6.1435 

1040  7 

-0-2 

Mo 

42 

95 

15.78 

5/2 

-0.9099 

277-4 

97 

9.60 

5/2 

-09290 

2833 

Pi 

46 

40r 

* 

22.23 

5/2 

-0.57 

174 

lr. 

49 

113 

4.16 

9/2 

5.4960 

931.0 

1  144 

Qi \§/Q  113=“ 

115 

9584 

9/2 

5.5072 

9329 

1  161 

=  10146 

Sb 

51 

121 

57.25 

5/2 

3.3417 

1019 

-053 

(/us'Ptii  = 

123 

42.75 

7/2 

2.5334 

551.8 

-0.68 

=  1.27475 

I 

53 

127 

100 

5/2 

2.7939 

8519 

-088 

Xe 

54 

131 

21.24 

3/2 

0.6868 

349 

-0.12 

Cs 

55 

133 

100 

7/2 

2.5642 

558.5 

-C.003 

Ba 

56 

135 

6.6 

3/2 

0.837 

425 

Q\nlQm  = 

137 

11.3 

3/2 

0.936 

476 

=  1.542 

La 

57 

139 

99.911 

7/2 

2.7615 

001-4 

0.5 

Pr 

59 

141 

100 

5  2 

3.8 

1130 

-0.054 

Nd 

60 

143 

12.2 

7/2 

-1.1 

220 

5=1.2 

145 

8-3 

7/2 

-0.69 

140 

sS1.2 

Sm 

62 

147 

15.07 

7/2 

-0.68 

147 

0.72 

149 

1385 

7/2 

-0.55 

119 

0.72 

Eu 

63 

151 

47.8 

5/2 

34 

1000 

—  1.2 

153 

52.2 

5/2 

1.5 

460 

—  2.5 

Yb 

70 

173 

16.08 

5/2 

—0.65 

198 

3.9 

Lu 

71 

175 

97.4 

7/2 

2.6 

570 

5.9 

Ta 

73 

181 

100 

7/2 

2.1 

460 

3.9 

Re 

75 

185 

37.1 

5/2 

3.1437 

958.6 

2.8 

187 

62.9 

5/2 

3.1760 

968  4 

2.6 

=  1.06 

Os 

76 

189 

161 

3/2 

9.6507 

3307 

08 

!  r 

77 

191 

38.5 

3/2 

0.16 

81 

1.5 

193 

61.5 

3/2 

0.17 

86 

1.5 

Au 

79 

197 

100 

3/2 

0.136 

69.1 

0.56 

Hg 

80 

201 

13.24 

3/2 

-0-607 

308 

05 

Bi 

83 

209 

100 

9/2 

4  0389 

684.3 

-0-4 

"T 
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2.2  GEOCHEMICAL  CONSIDERATIONS 


A  substantial  proportion  of  natural  nuclides  exhibit  NQR.  These  are  listed 
in  table  2,  which  was  taken  from  the  recent  comprehensive  book  by  Biryukov, 
et  al  (1969).  Most  of  these  are  relatively  rare  in  the  earth's  crust,  but 
a  few  are  among  the  petrogenic.  or  "lithophilic"  elements  found  in  quantity 
in  igneous  rocks  (Rankama  and  Sahama,  1950)  as  summarized  in  table  5. 


Lithophilic 

element 

Table  3.  The  most 

abundant 

NQR-active  isotopes 

Isotopic 
abundance  in 
igntous  rocks 

Percent  in 
igneous  rocks 

NQR 

isotope 

Abundance , 
percent 

Aluminum 

8.13% 

27A1 

100% 

8  13% 

Sodium 

2.8% 

23Na 

100% 

2.8% 

Potassium 

2.59% 

39K 

93.4% 

2.43% 

Magnesium 

2.09% 

Z5Mg 

11.5% 

0.24% 

Potassium 

2.59% 

41k 

6.6% 

0.17% 

Manganese 

0.1% 

55Mn 

100% 

0.1% 

Oxygen 

46.6% 

“  -17o_  ‘ 

0.041% 

0.019% 

Calcium 

3.63% 

43ca 

0.15% 

0.005% 

Others  below  0.1%  isotopic  abundance:  ^Ti,  35^]^  Up. 


The  nuclides  below  the  broken  line  in  table  3  may  not  be  as  high  in  isotopic 
abundance  (igneous  rocks)  as  some  of  the  minor  elements,  but  are  the  only 
important  NQR-active  constituents  of  major  sedimentary  rocks  such  as  sand¬ 
stone,  limestone  and  evaporites. 

It  is  unfortunate  that  the  most  abundant  rock-forming  elements,  oxygen  16  and 
silicon,  as  well  as  carbon,  are  not  NQR-active. 

The  most  abundant  minerals  in  igneous  rocks,  according  to  several  authorities 
(Randama  and  Sahama,  1950),  are  the  following: 

Feldspars  -  59.5%  to  67.5% 

Quartz  -  9%  to  l"7. 6% 

Mg,  Fe,  Ca  silicates  (hornblende, 

pyroxenes,  micas)  -  12  to  20% 

Ilmenite  -  1%  to  2% 

Apatite  -  0.4%  to  0.7%. 


Clearly,  the  aluminum  present  in  all  feldspars,  K,Na(AlSi30s)  or  Ca.BafA^S^Og)  , 
offers  the  best  possibility  as  an  NQR-active  element  widespread  in  natural 
rocks.  Secondary  possibilities  exist  with  the  sodium  and  potassiun  also 
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present  in  many  feldspars.  Fortunately,  nuclear  quadrupole  coupling  constants 
have  been  measured  in  two  feldspars,  albite  and  r.i. crocline,  as  a  by-product 
of  NMR  measurements  (Hafner  and  Hartmann,  1962).  These  are  listed  in  table  4, 
together  with  those  of  all  the  other  reasonably  abundant  minerals  for  which 
data  have  been  found  in  a  very  thorough  literature  search.  Unfortunately, 
most  NQR  work  in  minerals  has  been  done  on  gem-quality  semi-precious  stones, 
or  on  relatively  rare  minerals  of  arsenic,  bismuth,  or  copuer. 


2.3  OTHER  RESULTS  FROM  THE  LITERATURE  STUDY 

Most  of  the  basic  quadrupolar  parameters  for  minerals  listed  in  table  4  were 
obtained  from  experiments  other  than  pure  nuclear  quadrupole  resonance.  Hence, 
for  guidance  in  our  experiments,  it  has  been  necessary  to  calculate  most  of 
the  expected  resonance  frequencies  shown,  using  the  formulas  of  table  1. 

These  frequencies  provide  a  guide  for  laboratory  searches  for  the  resonances, 
although  discrepancies  of  up  to  ±  10  percent  are  expected  tc  be  found.  Some 
errors  arise  from  the  inexactness  of  the  formulas  of  table  1  where  n  is  large, 
others  from  variations  of  the  parameters  with  temperature,  (the  NQCC  generally 
decreases  at  higher  temperature),  and  still  others  from  variations  caused  by 
dispersed  impurities,  which  locally  distort  the  crystal  lattice  slightly 
(Schempp  and  Bray,  1970)  and  vary  among  different  samples. 

]t  will  be  noticed  in  table  4  that  there  are  at  least  two  resonances  for  2?A1 
in  each  mineral,  and  there  may  be  many  more  (Kyanite  has  8).  The  same  is  true 
of  17o.  The  feldspars  and  other  simple  alumino-si licates  have  all  their  27A1 
Dure  quadrupole  resonances  between  about  300  kilohertz  and  1  megahertz.  The 
170  resonances  between  1.47  and  5.36  MHz  for  ice  and  AI2O3  suggest  that  if 
the  oxygen  resonance  can  be  observed  in  quartz,  it  will  be  found  in  this  gen¬ 
eral  vicinity,  or  possibly  up  to  13  MHz  (Lucken,  1969).  No  published  NQR 
data  for  170  in  quartz  were  found,  despite  extensive  searching.  Sodium-23 
and  potassium-39  resonances  may  also  be  found  in  this  range  or  lower  (Lucken, 
1969),  judging  from  the  meager  data  available. 

These  frequencies  are  generally  considered  "too  low  for  convenient  study" 
(Schempp  and  Bray,  1970).  This  is  probably  because  the  absorption  line  widths 
and  separations  are  an  appreciable  fraction  of  the  center  frequency.  The  in¬ 
strument  most  favored  in  this  frequency  range  is  the  marginal  oscillator,  with 
frequency  modulation.  However,  "FM  has  not  been  successful  where  the  lines 
are  very  broad"  (Schempp  and  Bray,  1970).  Data  on  mineral  resonance  line 
widths  are  sparse  and  scattered.  Marfunin  (1965)  states  that  the  NQR  signal 
amplitude  is  proportional  to:  (resonance  frequency/line  width) ,  with  the 
usual  stationary  detection  methods.  This  means  that  the  detectability  of 
the  lines  depends  greatly  on  the  rjrity  and  crystal  perfection  of  the  specimen. 
For  this  reason,  most  NQR  studies  of  minerals  have  employed  optically  clear, 
semi-precious  stones.  Cloudy  specimens  have  not  given  NQR  lines.  However, 
Marfunin  (1965)  points  out  that  pulsating  techniques  are  not  subject  to  the 
line-width  limitation  on  amplitude,  and  that  by  their  use  NQR's  can  be  ob¬ 
served  in  crystals  where  they  would  be  otherwise  undetectable  (also  see  Smith, 
1971).  One  sophisticated  example  is  the  experiment  of  Slusher  and  Hahn  (1964), 
in  which  23sja,  35ci ,  and  37ci  pure  quadrupole  resonances  were  detected  by  a 
special  pulse  technique  in  salt,  NaCl,  despite  up  to  0.2  percent  of  lattice 
impurities  and  otherwise  cubic  symmetry. 
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Table  4.  Summary  of  selected  mineral  NQR  data 


NQCC 

Isotope 

_I__ 

Mineral 

(e2Qq/h),  MHz 

n 

Resonances,  MHz 

Reference 

2?A. 

5/2 

A 1  b  1 1  e 

3.29,  3.  37 

0.62,  0.634 

(0.655), (0.932) 

? 

M 

•  1 

Mir  rod  ine 

1.22 

0.21 

(0.506), (0.956) 

i 

II 

tl 

Spinel 

1.  72 

0 

(0.556), (1.12) 

j 

II 

II 

Topaz 

1.55,  1.67 

0.8,  0.38 

(0.334), (0.287) 
(0.434) ,(0.488) 

a,  b 

•  I 

•  I 

Corundum 

2.  393 

0 

0.359,  0.718 

c 

" 

•  1 

Spodumene 

2.950 

0.94 

0.751,  0.793 

d 

•  1 

II 

Beryl 

3.0Q3 

0 

0.464,  0.927 

d 

•  1 

II 

Chrysoberyl 

2.85,  2.846 

0.94,  0.7  6 

(0.63), (0.605) 
(0.81)  ,(0.80) 

e 

II 

Kyanlte 

3.7,  6.5 

0.89,  0.59 

(0.812)',  (1.27) 
(1.045)  ,(1.85) 

f 

II 

II 

II 

9.4,  10.0 

0.38,  0.27 

(1.43), (1.62) 

(2. 74), (2. 96) 

f 

II 

II 

Andal usi to 

5.9,  15.6 

0.69,  0.08 

(1.22),  (2.32) 
(1.66), (4.65) 

f 

" 

II 

A  (metal) 

••0.6  ) 

- 

<0.94,  <0.188 

h 

170 

•  • 

HjO  ice) 

9.8),  6.66 

0.4,  0.93 

(1.71), (1.47) 

(2. 86), (1.89) 

1,  g  >  k 

" 

II 

Corundum 

18.9 

0.62 

(3. 76), (5. 36) 

m 

2hta 

1/2 

Alblte 

2.64,  2.62 

0.25 

(1.29) 

9 

23Nn 

" 

NaCl  (gas) 

5.40 

- 

(2.70) 

8 

39K 

II 

Kcl  (gas) 

5.67 

- 

(2.83) 

g 

35C1 

A1C1  (gas) 

8.8 

- 

(4.4) 

g 

2)1 

1 

1U0  (gas, ice) 

0.318,  0.214 

0.06,  0.1 

(0.234;, (0.155) 
(0.243), (0.166) 

k 

7 1.1 

3/2 

Spodumene 

0.0757 

0.79 

(0.0417) 

d 

NOTES : 

1.  Some  crystals  have  the  active  Isotope  at  more  than  one  type,  of  site 
in  the  unit  cell,  hence  have  several  values  of  r  and  NQCC. 

2.  Underlined  values  are  alternatives,  given  in  secondary  references. 

3.  Parentheses  denote  my  calculation,  from  data  given. 
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Pen'kov  and  Sofin  (1967)  say  that  the  intensity  of  an  NQR  signal  is 
proportional  to  (gyromagnetic  ratio)2  for  the  active  nucleus.  As  can  be  seen 
in  table  2,  27A1,  23Na  and  7Li  have  large  gyromagnetic  ratios,  and  should 

give  NQR  signals  7  to  16  times  as  large  as  -Cl,  which  is  one  of  the  most 

popular  and  easily  detected  isotopes.  This  is  encouraging.  The  signals  from 
170,  39k,  and  2H  should  be  more  nearly  comparable  to  that  from  33C1. 

Pressure  applied  to  a  crystal  generally  increases  the  NQCC  and  the  resonant 
frequencies  (Schempp  and  Bray,  1970).  It  may  also  reduce  the  signal  amplitude 
The  magnitude  of  the  pressure  coefficient  appears  to  depend  upon  the  compres¬ 
sibility,  the  degree  of  covalent  bonding,  and  the  volume  dependence  of  thermal 
vibrations  in  the  crystal.  The  papers  of  Kushida  et  al ,  Armstrong  and  Jeffrey, 
Early  et  al,  and  D'Alessio  and  Scott  listed  in  part  II-C  of  appendix  2  discuss 
the  details.  Near  room  temperature,  some  typical  pressure  coefficients  of 
frequency  3v/3P  are: 

for  33C1  in  KCIO3,  +0.46  Hz/psi  for  v  =  23  MHz. 

for  63Cu  in  CugO,  +  2.47  Hz/psi  for  v  =  26  MHz. 

for  2^Na  in  NaN'03,  -  0.33  Hz/psi  for  v  =  0.16  MHz. 

In  the  last-named,  the  coefficient  becomes  zero  at  about  200°C  and  positive 
above  that  temperature. 

There  is  a  Zeeman  splitting  and  weakening  effect  on  NQR  lines  by  external 
magnetic  fields.  This  effect  can  be  used  to  modulate,  and  thus  identify,  NQR 
absorptions,  with  advantages  in  stability  over  frequency  modulation,  the  tech¬ 
nique  generally  used  in  steady-state  NQR  work  (Smith,  1971). 

Powdered  or  polycrystalline  samples  appear  to  be  generally  satisfactory  for 
NQR,  except  in  ductile  (metallic)  materials  in  which  cold-working  residual 
strains  may  be  created.  However,  "surface  defects"  may  cause  weak  spurious 
signals  in  finely  powdered  samples  (Slusher  and  Hahn,  1964) .  Annealing  should 
be  done  in  vacuum  to  prevent  absorption  of  gaseous  impurities  from  the  atmos¬ 
phere.  Thermal  annealing  at  temperatures  approaching  melting  can  improve  the 
NQR  signal  in  some  materials;  in  others  a  partial  phase  transition  apparently 
occurs  which  introduces  much  disorder  in  the  crystal  an£  destroys  the  signal 
(Brun,  Hafner  and  Hartmann,  1960).  Hence,  de-stressing  of  minerals  by  crush¬ 
ing  appears  to  be  preferable  to  heating,  if  it  should  be  required. 


2.4  FIRST  NQR  EXPERIMENTS 

The  experiments  began  with  a  prac+ical  search  for  pure  nuclear  quadrupole 
resonances  in  common  rock-forming  minerals.  A  marginal-osci llator-type  cw 
spectrometer  built  in  1969  at  Teledyne  Geotech  for  NQR  studies  of  nitrogen 
compounds  (Cook,  1970)  was  restored  to  operation  and  several  improvements 
were  added.  The  oscillator  and  buffer  amplifier  circuit  used,  taken  from 
Gill  et  al  (1967),  is  shown  in  figure  2.  It  was  followed  by  a  conventional 
audio  band-pass  amplifier  (Gill  et  al,  1967,  figure  2)  passing  about  L.00  to 
1500  Hz.  The  output  was  displayed  on  an  oscilloscope,  whose  sweep  was  line- 
synchronized.  The  Varicap  frequency  modulator  was  driven  at  60  Hz  by  a  6.3  v 
filament  transformer.  A  50  MHz  frequency  counter  driven  by  the  E180F  rf 
buffer  was  used  to  measure  frequencies.  Oscillator  harmonics  and  waveforms 
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were  monitored  with  a  communications  receiver  and  a  50  MHz  oscilloscope. 
Interchangeable  rf  tank  coils  were  provided  having  from  1  turn  to  60  turns  of 
2-1/4  inches  inside  diameter,  to  tune  from  1.5  to  58  MHz  with  the  100-100  pfd 
twin  variable  tuning  condensers.  Unfortunately,  this  circuit  would  not  oscil¬ 
late  from  about  4  to  10  MHz,  despite  several  kinds  of  feedback  adjustments. 
Oscillation  above  this  gap  was  of  low  amplitude  and  apparently  occurred  in  a 
different  mode. 

The  equipment  was  checked,  as  in  1969,  with  a  sample  of  hexamethy lene  tetra- 
mine  (Q^jg  N4 ,  whose  nitrogen-14  resonances  near  5.3  MHz  were  easily  detected. 
In  pure  samples,  only  one  resonance  is  supposed  to  occur.  Our  commercial- 
grade  (Eastman)  sanple,  however,  probably  contains  homologous,  slightly  dif¬ 
ferent  compounds,  since  we  have  always  observed  more  than  a  dozen  closely- 
spaced  absorption  "lines."  The  spectrometer  was  found  to  be  working  and 
adjustments  were  made  for  strongest  signals. 

Samples  of  Ottowa  sand  were  obtained  from  a  local  cement  manufacturer.  This 
sand,  used  as  a  laboratory  standard,  is  quarried  in  Ottowa,  Illinois.  Instead 
of  the  resonance  at  45.25  MHz  in  this  sand  reported  by  0.  K.  Kowallis  (1967), 
we  observed  a  group  of  absorption  lines  just  belcw  22  NUT z ,  and  none  near  45 
MHz.  This  may  indicate  that  the  resonance  found  by  Kowallis  was  really  at 
22.62  MHz,  and  that  the  frequency  was  misread  in  his  experiment  because  of  a 
prominent  second  harmonic  from  the  oscillator.  Attempts  were  made  to  find 
these  22  MHz  resonances  in  large  crystalline  quartz  specimens  from  Brazil  and 
from  Arkansas,  without  success.  This  result,  together  with  the  likelihood 
that  170  nuclear  quadrupole  resonances  in  other  simple  compounds  occur  below 
1  MHz  (table  4),  indicates  that  the  resonance  in  the  Ottowa  sand  is  an  EPR  or 
EFR  impurity  resonance,  and  does  not  arise  from  the  oxygen- 17  component  of 
the  quartz. 

An  attempt  was  made  to  remove  the  impurity-carrying  grains  from  the  Ottowa 
sand  by  "combing"  it  with  a  strong  magnet.  Some  grains  were  thus  removed, 
but  the  remaining  sand  still  gave  the  resonances,  and  the  separated  fraction 
did  so  only  in  the  same  proportion  to  its  volume  as  the  original  sample.  It 
was  concluded  that  the  resonances  were  not  caused  by  separate,  palpable  grains 
of  magnetite  or  other  ferrimagnetic  minerals.  A  clue  to  their  origin  is  that 
of  the  two  grades  of  Ottowa  sand  obtained,  the  coarse  #190  gave  much  stronger 
signals  than  the  fine  #109. 

Empirical  manual  searches  for  resonances  were  made  between  1.5  and  4  MHz  on 
1  lb  samples,  which  nearly  filled  the  coils,  of  clear  quartz,  clear  calcite, 
opaque  ruby  (corundum),  opaque  orthoclase  feldspar,  and  red  granite.  No 
strong  resonances  were  found,  except  spurious  ones  from  the  equipment  or  en¬ 
vironment  (radio  stations). 

Upon  completion  of  the  literature  study,  it  became  clear  that  most  of  the 
pure  quadrupole  resonances  sought  should  be  found  below  1  MHz  (table  4)  ;  also, 
they  probably  require  a  very  narrow  equipment  bandwidth,  as  provided  by  very 
slow  automatic  frequency  scanning  and  phase  detection,  to  be  separated  from 
noise.  Hence,  the  equipment  was  modified  to  operate  at  frequencies  down  to 
0.31  MHz,  was  fitted  with  a  slow  motor-drive  for  the  tuning  condensers  (10  min¬ 
utes  or  5  hours  per  scan)  and  a  limit  switch.  Phase  detection,  low-pass  fil¬ 
tering  and  strip-chart  recording  were  also  added.  A  block  diagram  of  the  final 
spectrometer  is  given  in  figure  3.  Systematic  searches  of  the  region  between 
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0  44  and  1.2  MHz  were  made  on  1  lb  specimens  of  clear  quartz,  Madagascar  ruby, 
synthetic  ruby  and  sapphire  (all  corundum),  calcite,  and  orthoclase  feldspar 
without  result.  In  the  course  of  these  searches,  many  combinations  of  fm 
modulation  excursion,  scan  rate,  field  strength,  oscillator  stability  setting, 
etc.,  were  made.  The  second  harmonic  of  the  modulation  frequency  (60  Hz  or 
400  Hz)  was  generally  used  for  the  best  compromise  between  sensitivity  and 
faise  amplitude  modulation.  A  small  tuned  LC  circuit  could  be  loosely  coupled 
to  vhe  specimen  coil  in  place  of  the  mineral  sample  for  checking  the  equipment 
response;  its  absorption  peak  was  much  broader  than  the  lines  seen  with  hexa- 
methylene  tetramine,  however. 

Careful  cearches  were  made  especially  for  the  two  known  resonances  of  27A1  in 
corundum  at  0.359  and  0.718  MHz,  as  revealed  by  the  splitting  of  the  NMR  line 
of  27ai  at  high  magnetic  fields  into  five  components  by  the  quadrupole  effect. 
Nothing  was  found  except  spurious  apparatus  resonances  which  persisted  with 
the  sample  removed.  (Some  of  these  are  attributed  to  the  piezoelectric  mica 
in  the  vacuum  tubes.) 

2.5  NQR  CONSULTATIONS  AND  ADVANCED  EXPERIMENTS 

Our  lack  of  initial  success  in  locating  suitable  nuclear  quadrupole  resonances 
in  samples  of  com  ..on  rock-forming  minerals  using  a  regenerative  C.W.  spectrome¬ 
ter  led  us  to  seek  the  advice  of  full-time  specialists  in  NQR.  Dr.  Thomas  A. 
Scott  of  the  University  of  Florida,  who  wrote  summaries  of  worldwide  NQR  work 
in  1970  and  1971  for  the  Magnetic  Resonance  Review,  was  contacted  first. 

Through  him,  other  specialists  were  located.  Dr.  Scott's  C.W.  spectrometers, 
used  in  studying  the  structures  of  nitrogen  compounds,  were  found  to  be  highly- 
developed  devices  which  had  been  refined  over  several  years.  Such  refinement 
was  of  course  impossible  for  us  to  achieve  in  the  present  short  program.  Dr. 
Scott  knew  of  no  NQR  work  on  27ai  or  1?0,  other  than  the  references  already 
located  (see  bibliography),  but  he  gave  much  general  advice  on  NQR  technique 
and  the  state  of  the  art. 

In  a  visit  to  Dr.  Wm.  Moulton  at  Florida  State  University,  we  were  advised 
that  ] ulse  spectrometers  are  far  more  sensitive  than  C.W.  spectrometers,  and 
were  shown  a  demonstration  of  the  *'C1  resonance  in  KCLO3  on  a  recently- 
designed  pulse  spectrometer.  Dr.  Moulton  advised  us  that  NQR  was  a  much  less 
promising  method  of  measuring  pressure  effects  in  minerals  than  are  EPR  (of 
colored  minerals)  or  EFR  (in  trace  minerals  of  "magnetically  ordered"  char¬ 
acter  such  as  magnetite).  He  and  his  staff  emphasized  that  the  scarcity  of 
170,  which  is  only  0.037%  of  all  oxygen,  the  absence  of  any  NQR  effect  in  the 
really  abundant  lithophilic  elements,  and  the  complexity,  variable  composition 
and  opaque  character  of  the  feldspars  and  other  silicates  containing  substan¬ 
tial  amounts  of  27A1,  should  discourage  any  attempt  to  employ  NQR  methods. 
Judging  from  subsequent  experience,  this  advice  was  correct. 

2.5.1  Pulse-Method  Searches 

Despite  Dr.  Moulton's  advice,  we  decided  to  explore  fully  the  quantitative 
limitations  of  NQR  in  rocks,  in  accordance  with  contract  provisions,  because 
of  the  broad  applicability  of  a  method  using  ^70  or  27A1,  and  because  of  the 
convenience  in  not  having  to  provide  an  aligning  magnetic  field.  Arrangements 
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were  made  with  Dr.  Slade  Carr  of  the  Florida  A  f,  M  University  to  search  for 
resonances  in  several  crushed  mineral  specimens  supplied  by  us,  using  an  ad¬ 
vanced  pulse  spectrometer.  Dr.  Carr  had  a  reputation  as  an  expert  searcher 
with  this  instrument.  Unfortunately,  it  was  found  incapable  of  being  tuned 
to  frequencies  below  2  MHz.  Empirical  searches  were  performed  from  2  to  4  MHz 
on  a  few  specimens,  and  from  4  to  8  MHz  on  specimens  of  13  minerals,  including 
synthetic  sapphire  and  ruby  (corundum),  natural  grey  corundum  and  ruby,  quartz, 
calcite,  orthoclase  and  albite  feldspars.  Nothing  was  found. 

Arrangements  were  made  with  Dr.  Tonis  Oja  at  the  University  of  Denver  to  per¬ 
form  searches  on  his  pulse  spectrometer  of  recent  design  (figure  4)  (Oja,  1972) 
which  was  capable  of  operating  at  frequencies  from  0.50  to  5.4  MHz.  Five  days 
of  tedious  searching  by  the  Principal  Investigator  and  by  Dr.  Oja  yielded  no 
2^A1  or  1'0  resonances  from  carefully  prepared,  crushed  50-gram  samples  of 
various  corundums  and  feldspars  purchased  from  Linde  and  from  Filer's  Mineral 
Supply  Co. 

The  searches  were  concentrated  around  the  re.sonant  frequencies  given  in  table  4, 
at  both  room  temperature  and  in  liquid  nitrogen  (7"7°K) .  The  two-pulse  method 
was  used,  employing  a  "90  degree"  pulse  and  a  "ISO  degree"  pulse  as  shown  in 
figure  5.  The  areas  of  these  pulser  were  carefully  adjusted  so  that  the  sec¬ 
ond  pulse  was  just  twice  the  first,  which  was  determined  by: 


t  H  =  -  /y  V (I  jm|)  (I+|m|+l) 


(3) 


The  pulse  length  t  and  field  strength  H  thus  depend  upon  the  gyromagnetic 
ratio  y  and  spin  I  of  the  nucleus  being  sought  (see  table  2) .  The  quantum 
transition  number  m  is  1/2  for  the  lower  frequency  and  is  3/2  for  the  upper 
frequency.  Since  the  r.t.  magnetic  field  strength  H  could  not  be  measured 
with  the  equipment  available,  it  was  estimated  from  the  resonance  voltage  E 
across  the  sample  coil,  which  was  part  of  a  parallel  resonant  ciicuit  (figure  4), 
according  to  the  following  formulas: 


i  =  E  /yj  X2  +  R2  =  EwCQ  /^Q2  +  1  ^  EwC  (4) 


H  =  Ni/10  L  -  KEwC/lOL  (5) 


Here  N  and  L  are  the  total  turns  and  length  of  the  coil,  C  the  tuning  capaci¬ 
tance  at  angular  frequency  <*>,  as  measured  on  a  Q-meter  for  each  frequency  and 
each  coil  before  installation,  and  E  is  the  r.f.  voltage  (peak  to  peak)  as 
measured  across  the  coil  with  an  oscilloscope.  These  formulas  permitted  the 
correct  adjustments  of  the  product  tE  to  be  made  for  each  of  the  several  coils 
required  for  the  various  resonance  ranges  searched.  The  coil  Q's  were  all 
above  100,  so  that  the  approximation  in  formula  (4)  is  permissible.  The  cal¬ 
culated  product  tE  was  checked  against  optimum  values  found  empirically  for 
nitrogen  and  niobium  compounds  which  gave  good  signals  in  the  equipment  (with 
proper  compensation  for  the  I,  y  and  m's  of  these  elements).  It  was  only 
necessary  to  make  tE  correct  to  within  about  a  factor  of  2. 
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A  sealed  specimen  of  aluminum  tribromide  was  also  tried  in  the  pulse  spectrome¬ 
ter.  This  specimen  had  given  a  strong  ^A1  NQR  signal  in  a  CW  spectrometer 
some  years  before.  No  resonances  were  found  by  the  pulse  method,  ihe  equip¬ 
ment  was  rechecked  on  some  nitrogen  compounds,  with  good  results;  in  fact,  two 
new  resonances  in  an  organic  compound  never  studied  before  were  discovered 
by  Dr.  Oja  after  about  an  hour's  searching.  Figure  6  illustrates  the  appear¬ 
ance  of  the  "decay  tail"  and  "echo"  signals  so  obtained,  which  were  sought  but 
never  seen  with  specimens  containing  aluminum. 

Because  of  the  possibility  of  sample  "saturation"  by  the  r.f.  field,  which 
could  rapidly  quench  the  signal,  the  searching  was  done  with  individual  pulse- 
pair  events,  applied  by  means  of  a  command  pushbutton  at  intervals  of  at  least 
1  second  (to  allow  time  for  thermalization  recovery  of  the  nuclear  spin  system). 
Between  events,  the  RF  oscillator  was  generally  changed  in  frequency  in  1  kHz 
steps.  About  each  10  kHz,  several  adjustments  on  the  transmitter  and  coupling 
network  had  to  be  retuned  to  preserve  the  tE  product  and  to  avoid  distortion 
of  the  pulses.  The  best  search  rate  achieved  was  only  100  kHz  per  hour.  In 
view  of  the  many  variables  involved,  all  the  possible  combinations  could  not 
be  thoroughly  explored  while  the  equipment  was  available. 

Our  failure  to  find  any  2</\l  resonances  was  attributed  to  one  or  more  of  the 
following  causes: 

a.  All  the  samples  may  have  been  unsuitable,  lacking  detectable  reso¬ 
nances  because  of  strains  or  impurities. 

b.  The  values  of  tE  used  and  the  frequency  ranges  searched  were  incor¬ 
rect,  despite  the  best  available  theoretical  guidance  and  much  care  in  calcu¬ 
lating  and  measuring  tE  by  alternative  independent  methods. 

c.  The  spin-lattice  relaxation  times  Ti  for  these  samples  were  much 
longer  than  1  second,  so  that  the  signals  did  not  appear  more  than  once  before 
saturation  occurred,  and  were  discounted  as  "noise." 

d.  The  values  of  T^  were  much  shorter  than  300  microseconds,  the  mini¬ 
mum  blanking  and  recovery  time  possible  for  the  particular  receiver,  hence, 
the  "decay  tail"  signal  was  over  before  the  receiver  was  working  and  there 
was  insufficient  nuclear  coherence  to  form  a  pulse  echo. 

To  test  hypothesis  (d) ,  arrangements  were  made  through  Dr.  Oja  to  use  a  new 
low-frequency  pulse  spectrometer  of  advanced  design  built  by  Mr.  Gary  Peterson, 
an  expert  specialist  with  MATEC,  Inc.  of  Warwick,  R.  I.  This  instrument  had 
a  receiver  "dead  time"  of  only  40  microseconds.  It  first  became  available 
during  the  January  1973  meeting  of  the  American  Physical  Society  in  New  York 
City,  where  it  was  being  demonstrated.  Twelve  hours  of  searching  for  the  A1 
resonances  in  various  corundums  and  feldspars  were  done  by  three  MATEC,  Inc. 
experts  and  by  the  Principal  Investigator.  For  much  of  this  work  a  digital 
transient-memory  and  display  unit  was  used  to  enhance  observation  of  the  short¬ 
lived  waveforms  displayed.  Both  room  temperature  and  liquid  nitrogen  tempera¬ 
ture  were  used.  Pulse  sequences  were  applied  as  much  as  30  seconds  apart  to 
avoid  possible  saturation  effects.  Nevertheless,  no  resonances  were  found. 
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2.5.2  Related  NMR  Findings 


In  the  meantime,  a  publication  was  obtained  from  Switzerland  (Hartmann,  1963) 
giving  details  of  NMR  experiments  on  natural  feldspars .showing  that  many  such 
minerals  can  show  the  quadrupole  splitting  effect  in  ^7A1  quite  clearly,  as 
seen  in  figure  7.  After  a  lengthy  series  of  inquiries,  arrangements  were  made 
with  Dr.  Irwin  Weisman  at  the  National  Bureau  of  Standards  to  run  NMR  spectra 
of  samples  furnished  by  us,  to  determine  if  they  exhibited  quadrupole  splitting 
and  thus  to  test  hypothesis  (a)  above.  He  immediately  found  the  quadrupole 
splitting  in  the  Linde  synthetic  sapphire  (figure  S) ,  but  only  weak  hints  of 
splitting  in  two  natural  feldspars.  A  major  reason  was  probably  that  in  the 
crushed  samples  furnished,  the  random  crystal  orientations  smeared  the  posi¬ 
tions  of  the  four  satellite  lines  too  much,  since  their  separation  from  the 
central  peak  (figure  7)  varies  greatly  with  crystal  orientation  in  the  NMR 
method  (Hartmann,  1963)  (crystal  orientation  affects  only  the  amplitude  of 
pure  NQR  signals,  hence  crushed  samples  give  a  signal  about  half  as  large  as 
an  optimally-oriented  single  crystal). 


Arrangements  could  not  be  made  to  have  Dr.  Weisman  's  initial  results  extended 
at  the  Bureau  of  Standards,  because  of  administrative  difficulties.  Attempts 
were  made  to  test  single  crystals  and  additional  samples  in  a  standard  high- 
resolution  NMR  spectrometer  at  an  institution  in  Texas,  without  success;  its 
bandwidth  and  modulation  were  far  too  narrow  for  the  wide  spectral  lines  found 
in  solids.  Calls  to  universities  and  major  government  and  industrial  labora¬ 
tories  throughout  the  nation  failed  to  locate  another  wide-line  NMR  spectrome¬ 
ter  equipped  to  operate  in  the  27^1  region  (which  is  near  that  of  carbon- 1 3) . 
However,  it  was  concluded  that  hypothesis  (a)  above  had  been  essentially  dis¬ 
proved  the  synthetic  sapphire,  at  least,  should  give  a  pure  NQR  signal. 

2.5.3  Final  C.W.  Experiments 

Finally,  arrangements  were  made  with  Dr.  Thomas  A.  Scott  at  the  University  of 
Florida,  for  his  staff  to  perform  pure  NQR  searches  on  some  mineral  specimens, 
using  his  advanced  CW  equipment  with  the  summation  of  repeated  frequency  scans. 
We  reasoned  that  whatever  the  errors  in  adjustment  (hypothesis  b)  or  peculiari¬ 
ties  of  27ai  (hypotheses  c  and  d)  which  prevented  its  detection  in  the  pulse 
experiments,  these  should  not  affect  the  CW  method. 

After  initial  difficulties  from  piezoelectric  behavior  of  the  synthetic  sap¬ 
phire,  zeeman  (magnetic  field)  modulation  was  used.  Radii.  station  interfer¬ 
ence  was  eliminated  by  careful  shielding.  Weak  resonances  were  finally  found 
at  0.363  MHz  and  0.722  MHz,  within  1%  of  the  predicted  values.  Figure  9  shows 
the  accumulated  result  of  128  passes  through  the  lower  resonance.  A  compari¬ 
son  recording  with  the  resonance  destroyed  by  an  inhomogeneous  d.c.  magnetic 
field  is  also  shown.  The  resonance  width  is  about  7  kHz  and  the  signal-to- 
r.oise  amplitude  ratio  is  about  2.  All  tests  were  made  at  room  temperature. 

Samples  of  an  albite  from  Gunnison,  Colorado,  and  a  miorec.line  from  Sonora, 
California,  also  gave  repeatable  but  very  weak  signals  it  frequencies  within 
3%  of  the  values  predicted  by  table  4.  Table  5  below  summarizes  all  results. 
These  results  required  a  week's  time  to  acquire,  in  a  well-equipped  laboratory 
specializing  in  low-frequency  NQR,  by  skilled  experts  with  years  of  experience 
(in  particular,  Paul  Canapa) .  They  are  surely  representative  of  the  best  that 
can  be  done  at  this  time. 
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(a)  Albite  single  crystal, 

from  Amelia  Court  House,  Va. 


(b)  Twinned  polycrystalline 
niicrocline  from  Crystal 
Peak,  Colorado 


(c)  Microcline  single  crystal 
from  Hull,  Quebec. 


w 


9,0  AS  8,0  7S  1C 


Figure  7.  NMR  spectra  of  three  feldspars  at  9  MHz,  with 

lines  from  27^1  and  23\'a  (shaded) ;  from  Hartmann 
(1963).  Optimum  orientations. 
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363  kHz 


Figure  9.  Averaged  recording  from  128  scans  of  pure  NQR  signal  from  synthetic  sapphire 


Table  5.  C.W.  NOR  Results 


No.  of  Frequency,  Line  width.  Signal/ 


Mineral 

scans  summed 

MHz 

MHz 

noise 

Synthetic 

j  128 

0.363 

0.007 

2 

sapphire 

(  100? 

0.722 

0.010 

3? 

Albite 

64 

0.955 

0.010 

1  •  5  (  ?) 

Microcline 

156 

0.925 

0.018 

1 • 5  (?) 

Throughout  the  early  empirical  searches  with  CW  and  pulse  spectrometers,  large 
single-crystal  specimens  (up  to  1  lb)  of  clear  calcite,  quartz,  and  corundum 
(Linde  synthetic  sapphire  and  ruby  scrap)  were  tested  repeatedly  in  the  fre¬ 
quency  range  of  1  to  6  MHz,  where  pure  NQR  resonances  of  can  be  expected 
according  to  table  4.  Nothing  was  ever  found  which  depended  upon  the  presence 
of  the  sample  in  the  coil.  The  optimum  orientation  of  the  crystals  in  the 
coil  was  never  determined,  but  some  chunks  of  scrap  quartz  from  Brazil  and  the 
sapphire  fragments  were  shorter  than  the  coil  diameter,  so  were  tried  at  vari¬ 
ous  orientations.  The  clear  crystals,  large  size  and  high  oxygen  content  of 
the  samples  used  in  these  tests  should  have  aided  in  producing  a  detectable 
NQR  signal. 

During  the  pulsed-spectrometer  searches  at  the  University  of  Denver,  the 
nitrogen- 14  signal  from  a  30  ml  sample  of  CgHnN  was  measured;  a  60  mv  signal 
amplitude  was  obtained  from  about  2.3  x  1023  nuclei  of  N.  The  same_sample 
volume  of  quartz  was  estimated  to  contain  only  4.5  x  1020  nuclei  of  I/O,  be¬ 
cause  of  its  low  isotopic  abundance.  On  the  other  hand,  according  to  Pen'kov 
and  Sofin  (1967),  the  (gyromagnetic  ratio)^  for  17(3  is  3.5  that  of  ^N  (see 
table  2).  Since  the  ratio:  frequency/line  width  is  unknown  for  170,  its 

influence  on  signal  strength  cannot  be  estimated.  Consideration  only  of  the 
ratios  of  nuclear  abundance  and  the  square  of  the  gyromagnetic  ratio  led  to 
the  expectation  that  a  quartz  sample  should  have  produced  a  signal  amplitude 
of  only  0.4  millivolt  in  the  equipment  at  the  University  of  Denver.  Since  the 
noise  level  was  1  to  2  mV,  no  further  serious  searches  for  the  17(3  resonances 
were  made. 

Subsequently,  conversations  with  Prof.  P.  J.  Bray  of  Brown  University,  a  re¬ 
spected  senior  NQR  researcher,  indicated  that  the  370  resonances  had  never 
been  observed  except  in  "heavy  ice,"  obtainable  from  Israel  enriched  with 
100  times  or  more  the  0.037%  isotopic  abundance  of  0  found  in  all  natural 
materials.  Still  later,  a  reference  was  found  (Brun,  et  al,  1970)  to  what  was 
believed  to  be  the  first  observation  of  the  quadrupole  effect  in  170  without 
artificial  enrichment.  The  signal  was  enhanced  by  a  large  but  unknown  factor 
by  employing  the  Overhauser  effect  (Abragam  §  Borghini,  1964).  This  apparently 
can  be  done  only  in  a  standard  NMR  experiment  and  requires  in  addition,  a  very 
low  temperature  ii  the  liquid  helium  range  and  a  strong  transverse  r.f.  mag¬ 
netic  field  at  the  resonant  frequency  of  the  electron  system  (9  to  72  gigaherz) . 
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The  paramagnetic  chromium-ion  doping  in  the  ruby  specimen  used  was  apparently 
necessary  to  fulfill  the  requirement  of  an  energy-saturated  electron  system. 

In  the  experiment  of  Brun,  et  al  (1970),  new  values  of  the  quadrupole  coeffic¬ 
ients  of  i-n  corundum  were  found:  NQCC  =  2.167  MHz  and  n  =  0.517  at  1.9°K. 

These  should  be  slightly  different  at  room  temperature,  but  would  nevertheless 
indicate  that  the  pure  NQR  lines  should  be  found  at  0.406  and  0.622  MHz,  very 
far  from  the  values  given  in  table  4  on  the  basis  of  theoretical  calculations 
by  Hafner  and  Raymond  (1968).  It  is  not  surprising  that  we  never  saw  the 
upper  one  in  the  course  of  searching  the  frequency  range  between  0.5  and  1  MHz 
for  the  upper  27A1  resonance,  in  view  of  the  delicate,  complicated  and  arti¬ 
ficial  laboratory  procedure  required  in  the  experiment  of  Brun,  et  al  (1970). 
Neither  such  procedures,  nor  artificial  oxygen-isotope  enrichment  of  the  rock 
samples  are  possible  as  a  part  of  a  borehole  or  field  method  of  measuring  rock 
stress. 


2.6  CONCLUSIONS  REGARDING  THE  USE  OF  NQR 

The  outlook  is  now  rather  poor  for  the  eventual  use  of  the  NQR's  of  27A1,  2%a 
and  39k  in  feldspars  to  measure  rock  stress.  The  feldspar  NQR  signals  of  27A1 
appear  to  be  substantially  weaker  than  the  27A1  signal  of  pure  corundum  (AI2O3) , 
which  was  itself  barely  detectable  by  a  sophisticated  laboratory  technique. 

This  is  probably  because  of  the  lower  concentration  of  aluminum,  plus  colored 
impurities,  in  the  feldspars.  The  same  consideration  would  apply  to  23^a  and 

A  major  impediment  to  the  use  of  all  three  isotopes,  as  well  as  170,  appears 
to  be  the  comparatively  low  frequency  of  the  resonances;  most  radio-frequency 
spectrometric  methods  become  more  sensitive  at  higher  frequencies.  One  basic 
reason  for  this  is  that  the  voltage  induced  in  a  coil  by  nuclear  or  electronic 
magnets  is  proportional  to  the  frequency  of  thei.  emissions,  in  accordance 
with  Faraday's  law  of  electromagnetic  induction.  From  a  quantum  point  of  view, 
the  energy  of  interaction  per  particle  is  proportional  to  the  frequency. 

In  all  the  advanced  NQR  experiments  performed,  samples  of  50  gm  or  less  were 
used.  There  is  some  possibility  that  the  use  of  larger  samples  (e.g.,  5  kilo¬ 
grams)  would  enhance  the  NQR  signals  of  feldspars  enough  to  make  them  detec¬ 
table.  However,  a  practical  method  of  measuring  rock  stress  in  situ  will  in¬ 
volve  a  considerable  further  degradation  of  the  signal,  because  the  rock  sample 
is  not  all  feldspar,  and  it  cannot  be  inside  the  coil  but  only  adjacent  to  it 
or  surrounding  it.  Further  quantitative  experiments  would  be  needed  to  deter¬ 
mine  the  detectability  of  NQR  signals  from  large,  external  rock  samples  as 
found  in  a  drill  hole. 

The  outlook  for  using  l70  to  measure  rt ck  stress  is  substantially  worse  than 
that  for  2  Al  or  2^Na,  because  of  the  much  lower  abundance  of  this  isotope  in 
natural  rock-forming  minerals. 
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3.  NUCLEAR  MAGNETIC  RESONANCE  METHODS 


Conventional  NMR  has  received  by  fur  the  greatest  amount  of  research  of  any 
of  the  radio-frequency  spectroscopic  methods.  An  entire  research  journal, 
Nuclear  Magnetic  Resonance,  has  been  devoted  largely  to  this  subject  since 
1968.  The  nucleus  receiving  the  bulk  of  the  effort  has  been  the  proton,  be¬ 
cause  hydrogen  is  a  prominent  constituent  of  water  and  most  organic  compounds, 
and  because  it  gives  a  particularly  strong  NMR  signal.  Material  in  the  liquid 
state  is  best  for  NMR,  since  the  rapid  molecular  motion  in  liquids  averages 
out  the  magnetic  disturbances  from  neighboring  atoms  and  permits  sharp  reso¬ 
nances  of  the  nucleus  with  the  external  fields. 

NMR  occurs  only  in  the  presence  of  a  highly  uniform,  steady  external  magnetic 
"aligning"  field,  about  which  the  spinning  magnetic  nuclei  can  p^ecess.  The 
precession  couples  magnetically  with  coils  surrounding  the  sample,  whose  axes 
are  perpendicular  to  the  aligning  field,  as  shown  in  figure  10.  The  frequency 
of  resonance  is  given  by: 


f  =  uH0/hI 


(6) 


Here  11^ 

I 

u 


aligning  field  strength  in  gauss. 
Planck's  constant  =  6.62  x  10"27  erg  sec 
nuclear  spin  number 
nuclear  magnetic  moment 


For  example,  in  a  field  HQ  of  5000  gauss,  typical  of  those  generally  used, 
hydrogen  nuclei  resonate  at  21.3  MHz.  The  earth's  field,  however,  is  suf¬ 
ficiently  uniform  over  a  sample  to  be  used  for  alignment.  For  a  typical 
earth ' s-field  strength  of  0.5  gauss,  the  proton  resonance  occurs  at  2,130  Hz. 
Excellent  digital  magnetometers  exist  which  employ  NMR  in  this  way.  The 
sample,  of  highly  purified  water,  hydrocarbon,  or  alcohol,  must  have  a  volume 
greater  than  about  500  ml  to  obtain  sufficient  NMR  signal.  The  protons  are 
pre-aligned  by  a  strong  but  non-uniform  dc  field  from  a  pulsed  coil  before 
each  measurement  interval;  when  the  pre-alignment  field  is  removed,  the  pro¬ 
tons  remain  sufficiently  aligned  in  the  earth's  field  to  give  a  strong  NMR 
signal  for  about  1  second,  the  spin-spin  coupling  relaxation  time  T2.  This 
time  is  made  long  enough  by  having  a  very  pure  liquid  sample. 


3. 1  FINDINGS  OF  THE  NMR  LITERATURE  STUDY 

A  Swiss  group  has  demonstrated  the  feasibility  of  NMR  studies  of  several 
organic  liquids  in  the  earth's  field,  with  or  without  a  strong  pre-aligning 
field  (Bene,  1966) .  Through  successive  improvements  in  equipment,  ending  in 
that  shown  in  figure  11  (Merck,  et  al,  1967),  they  have  achieved  a  resolution 
of  0.06  Hz.  They  had  sufficient  sensitivity  to  detect  weak  satellites  of  the 
proton  resonance,  caused  by  interaction  with  rare  isotopes  such  as  carbon- 13, 
which  constitutes  about  1%  of  the  carbon  in  organic  liquids.  Figure  12  shows 
the  satellite  group  found  for  methanol,  CH3OH.  The  calculated  spectral  lines 
are  shown  below  as  vertical  bars. 
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Figure  10.  Schematic  diagram  of  Rollin' s  arrangement  for  observing 
nuclear  magnetic  resonance  (from  Andrew,  1954) 
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Figure  11.  Equipment  for  NMR  in  the  earth's  magnetic  field 
(from  Merck,  et  al,  1966) 
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Figure  12.  Proton  satellite  resonances  for  the  C  -  H  band 
in  CH^OH  (from  Bene,  1966) 
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Figure  13.  Resonance  line  shape  for  the  protons  in  polycrystalline 
gypsum  CaS04.2H20,  taking  into  account  the  interaction 
of  all  neighbor  atoms  (from  Andrew,  1954) 


The  frequencies  of  such  satellites  should  be  stress-sensitive,  sirce  they 
depend  upon  the  separation  of  the  protons  from  the  13C  atom  in  the  molecule. 

It  appears  that  a  similar  scheme  might  be  considered  for  NMR  measurements  of 
rock  stress  in  the  earth's  field.  Numerous  solids,  including  the  natural 
minerals  AI9O3  (corundum),  CaF2  (fluorite),  LiAl*(Si03)2  (spodumene),  BeA^Oq 
(chrysoberyl) ,  CaSC^ . 21^0  (gypsum),  H2O  (ice),  and  glass  have  been  investi¬ 
gated  structurally  (at  high  fields)  by  NMR  (Andrew,  1969).  Besides  hydrogen, 
about  10  nuclides  including  ^Al,  l?o,  29gj^  and  43^a  have  been  detected  in 
minerals  by  NMR.  Signal  strengths  are  low,  resonance  patterns  are  very  broad 
and  their  detailed  structure  is  badly  blurred,  particularly  for  polycrystal¬ 
line  samples,  as  shown  in  figure  13  for  gypsum.  Nevertheless,  it  has  been 
shown  that  the  mean-square  width  of  the  resonance  (its  "second  moment")  is 
often  proportional  to  r-3,  where  r  is  the  mean  separation  between  pairs  of 
resonating  nuclei,  hence,  the  NMR  pattern  width  is  very  sensitive  to  r,  which 
has  been  measured  to  a  precision  of  0.4  percent  via  NMR  pattern  widths.  Since 
r  changes  under  mechanical  stress,  rock  stress  can  be  measured  by  NMR,  in 
principle. 

In  many  solids  the  resonance  pattern  width  is  determined  largely  by  local 
magnetic  fields  from  neighboring  nuclei.  For  example,  the  gypsum  curve  of 
figure  13  consist. s  of  overlapping  resonances  from  the  two  protons  in  the  water 
of  hydration.  Fa  ;h  perturbs  the  magnetic  field  at  the  other  by  an  amount 
u/r3  =  ±  5  gauss,  equivalent  to  a  minimum  pattern  width  of  42.3  kHz,  which 
could  only  be  plotted  with  a  higher  center  frequency.  This  would  require  an 
external  aligning  field  of  perhaps  50  gauss,  obtainable  in  the  central  region 
of  a  4000  ampere-turn  loop  50  cm  in  radius,  for  example.  The  practicability 
of  such  a  scheme  requires  further  study,  since  unavoidable  nonuniformities  in 
the  aligning  field  might  broaden  and  weaken  the  gypsum  pattern  further. 

Line  widths  for  protons  in  32  minerals  are  tabulated  by  Marfunir.  (1965).  They 
range  from  0.1  gauss  (425  Hz)  in  kaolinite  to  11  gauss  (4.7  kHz)  in  bauxite. 
Most  would  require  a  strong  artificial  aligning  field  to  be  observed. 

In  part,  the  line  width  in  solids  appears  to  depend  in  a  complicated  way  upon 
the  spin-spin  and  spin- lattice  relaxation  time  T2  and  Ti,  which  in  turn  de¬ 
pend  on  composition,  particularly  on  traces  of  paramagnetic  impurities  (Bloem- 
bergen,  1961).  Generally  T2  is  short  and  Tj  long  compared  to  the  values  for 
liquids,  which  results  in  a  poor  signal-to-noise  ratio.  Andrew  (1969)  states 
that  the  NMR  spectrum  of  a  solid  can  rarely  be  presented  on  an  oscillograph, 
and  instead  its  first  derivative  is  usually  recorded  using  a  phase  sensitive 
amplifier. 

An  encouraging  precedent  for  earth's  field  NMR  in  solids  exists  in  that  the 
?Li  and  resonances  in  a  very  pure  sample  of  LiF  (with  a  vei  y  long  relaxa¬ 

tion  time)  have  been  measured  in  the  earth's  field  by  use  of  a  strong  pre¬ 
aligning  field  (Ramsey  §  Pound,  1951).  Although  the  pattern  was  about  200  kHz 
broad,  there  were  peaks  at  the  central  resonant  frequencies  of  these  elements. 


3.2  POSSIBILITIES  FOR  NMR  AS  A  METHOD  OF  MEASURING  ROCK  STRESS 

The  general  requirements  for  successful  NMR  of  natural  polycrystalline  minerals 
in  situ  would  appear  to  be  as  follows: 


1.  Pre-alignment  by  a  strong,  temporary  dc  field; 

2.  Large  sample  volume  (several  kilograms); 

3.  Narrow  resonance  pattern  width,  of  the  order  of  100  Hz. 

4.  Low  r.f.  power  level  to  avoid  sample  saturation. 

In  the  laboratory,  requirement  (3)  is  met  in  many  materials  by  working  at  very 
low  temperatures.  In  our  application  this  is  not  possible;  our  best  alterna¬ 
tive  is  to  select  an  "indicator"  mineral  with  inherently  sharp  resonances. 

We  concluded  that  rock-stress  measurement  by  NMR  may  be  difficult  to  develop 
and  experiments  would  not  be  attempted  until  the  possibilities  of  NQR  had  been 
explored.  Nevertheless,  NMR  offers  the  possibility  of  employing  several 
additional  minerals  as  stress  transducers,  since  there  are  many  elements  capa¬ 
ble  of  NMR  but  not  NQR,  such  as  hydrogen  and  silicon,  which  are  prominent  in 
natural  rocks.  The  possibility  of  using  NMR  should  be  left  open,  for  the 
present . 


3.3  ZERO-FIELD  NMR  OF  MAGNETICALLY  ORDERED  MINERALS  (ZNR) 

In  magnetically-ordered  (ferromagnetic,  ferrimagnetic,  or  antiferromagnetic) 
crystals,  there  is  a  strong  internal  magnetic  field  arising  from  the  electron 
spins  and  orbital  currents  coordinated  by  strong  coupling  between  nearby  atoms. 
In  this  field  any  nucleus  present  which  has  a  magnetic  moment  will  have  a 
built-in  resonance  frequency. 

Several  magnetically-ordered  minerals  occur  as  accessory  minerals  in  common 
rocks.  As  noted  on  page  8  of  this  report,  igneous  rocks  contain  an  average 
of  1  to  2%  of  ilmenite,  FeO-TiC^,  which  has  a  magnetic  susceptibility  of  0.03 
to  0.14  cgs  units  (Jakosky,  1950).  Magnetite,  Fe304,  has  a  susceptibility  of 
0.04  to  24  cgs  units  (Jakosky,  1950)  and  occurs  as  very  fine  particles  in  basalt 
and  other  dark  igneous  rocks.  It  provides  the  basis,  together  with  its  mix¬ 
tures  with  ilmenite  and  ulvospinel  (spinel  =  MgO»Al203),  of  the  science  of 
paleomagnetism  (Dunlop,  1973).  Magnetite  also  constitutes  from  20  ppm  to  0.2?o 
of  most  terrestrial  soils  (Cook  and  Carts,  1962),  so  presumably  occurs  also 
in  sedimentary  rocks  such  as  shales  and  dark  sandstones,  since  it  is  a  chemi¬ 
cally  resistant  mineral.  Other  "magnetic"  minerals  found  widely  distributed 
in  nature  include  maghemite,  (y  -  Fe203)suspected  in  basalts  (Dunlop,  1973)  and 
in  red  soils  and  rocks  (Cook  and  Carts,  1962),  and  some  of  the  garnets  (common 
black  garnet  =  3  CaO- Fe203’ 3Si02)  . 

The  joint  effect  of  these  minerals  in  producing  magnetic  susceptibility 
(Jakosky,  1950)  permits  the  estimation  of  the  equivalent  percentage  of  mag¬ 
netite  (assumed  susceptibility  =0.2  cgsu  from  Dunlop,  1973)  in  typical  major 
rock  types,  as  summarized  below: 
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Table  6.  Estimated  magnetite  content  of  rocks 


Rock 

Rock  susceptibility,  cgsu 

%  "magnetite" 

Basalt 

0.00068  to  0.0063 

0.3  to  3 

Granite 

0.00003  to  0.0027 

0.015  to  1.3 

Gneiss 

0.00001  to  0.002 

0.005  to  1. 

Schist 

0.000026  to  0.00024 

0.013  to  0.12 

Sh  ile 

0.00004  to  0.00005 

0.02  to  0.025 

Sandstone 

0.000005  to  0.000017 

0.0025  to  0.008 

Limestone 

0.000004  - 

0.002  - 

It  is  apparent  that  the  ZNR  of  magnetite  and  related  ferrimagnetic  minerals 
could  furnish  a  basis  for  a  widely  applicable  stress-measuring  method  for 
rocks,  providing  the  effect  is  strong  enough  and  has  a  suitable  pressure  co¬ 
efficient.  The  45.5  MHz  resonance  of  iron  particles  measured  by  Kowallis  in 
zero  applied  magnetic  field  was  a  ZNR  (Benedek,  1965).  That  particular  reso¬ 
nance  is  of  only  passing  interest  here,  since  free  iron  is  almost  never  found 
in  the  earth.  We  have  confined  the  small  amount  of  effort  on  ZNR  in  this 
program  to  magnetite. 

No  attempt  was  made  to  survey  the  extensive  published  literature  on  the  NMR 
of  magnetic  materials  (ZNR).  Most  of  it  deals  with  synthetic  ferrites  and 
garnets,  and  miscellaneous  inorganic  compounds  rare  or  absent  in  nature. 
Andrew  (1972)  lists  29  recent  references  to  this  work.  The  active  nuclei 
used  as  "probes"  within  the  various  compounds  included  hydrogen,  vanadium, 
manganese,  fluorine,  iron-57,  europium,  etc.  Very  little  of  this  work  is  of 
any  interest  to  this  program. 

From  a  reference  (Portis  and  Lindquist,  1965)  given  us  by  Dr.  Moulton  at  the 
Florida  State  University,  it  was  learned  that  zero-field  resonances  had  been 
found  in  magnetite  at  67.621  and  63.549  MHz,  and  haematite,  a  -  Fe2C>3,  reso¬ 
nates  at  71.5  MHz  at  room  temperature.  The  signal-to-noise  ratio  of  the 
upper  magnetite  resonance  was  reported  to  be  "very  good"  and  the  line  "sharp" 
(Ogawa  and  Morimoto,  1962);  the  temperature  coefficient  of  the  resonances  of 
both  magnetite  and  natural  haematite  was  found  to  be  3.6  x  10'Vdeg  C,  with 
frequency  decreasing  as  temperature  is  raised.  Through  inquiries,  we  learned 
of  recent  precise  work  on  magnetite  at  the  Naval  Re-earch  Laboratory  and  ob¬ 
tained  a  preprint  of  a  new  paper  giving  the  resonance  line  shapes  (Rubinstein, 
Stauss  and  Bruni,  1973).  These  are  reproduced  for  room  temperature  in  fig¬ 
ure  14.  At  liquid  helium  temperatures,  the  upper  line  narrows  and  can  be 
resolved  into  four  components,  each  only  about  0.3  MHz  wide  (that  is,  0.4%  of 
the  frequency).  The  room-temperature  line  is  about  1.2  MHz  wide,  but  the 
position  of  its  peak  might  be  measured  to  within  0.25  MHz  or  0.35%  of  the 
frequency.  However,  George  Stauss  of  NRL  explained  that  the  line  shapes 
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ROOM  TEMPERATURE 
SPIN  ECHO  SPECTRUM 
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Figure  14. 


Figure  15. 


Room-temperature,  zero-field,  pulse-echo  NMR  spectrum  of 
57Fe  in  magnetite  (from  Rubinstein,  Stauss  and  Bruni,  1973) 


CC  2  4  6  8  10  \2 

H0  CnOe) 


Antiferromagnetic  resonance  frequency  versus  applied 
field  H0  for  aFe203,  (from  Foner,  1963) 
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shown  in  figure  14  were  smeared  by  the  pulse  width  of  5  microseconds  used, 
and  that  a  c.w.  spectrogram  should  give  better  resolution. 

3.3.1  Experiments 

The  marginal-oscillator  c.w.  spectrometer  at  Teledyne  Geotech  was  again  modi¬ 
fied  to  operate  at  frequencies  from  47.5  to  72  MHz,  and  samples  of  powdered 
magnetite  and  ilmenite,  diluted  with  50%  or  more  of  sugar  or  salt  (to  insure 
good  penetration  of  the  r.f.  magnetic  field)  were  empirically  tested  for  reso¬ 
nances.  Nothing  was  found  but  spurious  equipment  resonances  which  could  not 
be  eradicated.  Oscillations  were  very  weak,  so  the  frequency  counter  and 
oscilloscope  were  unreliable.  Since  these  tests  were  made  near  the  end  of 
the  program  period,  no  attempt  was  made  to  build  a  new  spectrometer  more  suita¬ 
ble  for  this  frequency  range,  and  to  insure  that  the  r.f.  field  was  below  the 
rather  low  saturation  threshold  for  magnetite. 

3.3.2  Further  Consultations 

In  conversations  with  several  research  people  familiar  with  ZNR  in  magnetite, 
it  developed  that  they  considered  these  resonances  weak  and  difficult  to  de¬ 
tect;  the  chief  reason  is  that  ^7Fe  comprises  only  2.45%  of  the  iron  present. 

No  published  measurement  of  the  pressure  coefficient  was  known.  George  Stauss 
attempted,  at  our  request,  to  assemble  equipment  to  measure  it,  but  was  unable 
to  do  so  in  the  short  time  available,  especially  since  magnetite  work  was  be¬ 
ing  terminated  at  NRL.  However,  it  was  learned  that  the  pressure  coefficients 
of  other  "hyperfine  resonances"  amount  to  0.01  to  0.1%  of  the  frequency  per 
ki lobar.  (The  pressure  coefficient  for  haematite  may  have  been  published  by 
David  Anderson  in  the  Physical  Review  about  1966,  but  was  not  looked  up.) 
Finally,  preliminary  data  were  obtained  on  the  pressure  effect  on  three  in¬ 
ternal  magnetic  fields  of  magnetite  at  room  temperature,  as  obtain  eel'  by  Dr. 
Nicholas  Halasa  at  the  University  of  Illinois  by  Mossbauer  Methods  (private 
communication) .  These  field  changes,  which  imply  a  decrease  in  the  upper 
(tetrahedral)  resonance  frequency  with  increasing  pressure,  average  approxi¬ 
mately  2.2%  for  100  kilobars  (1,400,000  psi) .  This  is  a  disappointingly  small 
effect,  and  is  within  the  range  stated  above  for  hyperfine  resonances. 

3.3.3  Conclusions  on  ZNR 

If  it  can  be  assumed  that  under  ideal  conditions,  the  frequency  of  the  upper 
magnetite  resonance  can  be  measured  to  a  precision  of  0.08  MHz  or  0.12%,  this 
would  mean  a  pressure  resolution  of  5.5  kilobars  or  76,000  psi.  This  would 
not  be  useful  for  practical  rock  stress  measurements.  It  would  be  necessary 
to  postulate  measurements  to  a  precision  of  1  kilohertz  or  better  to  resolve 
1,000  psi,  which  would  be  an  excessive  and  indefensible  extrapolation  of 
present  technology.  In  view  of  the  very  low  concentrations  of  magnetite  in 
most  rocks  (table  6),  the  low  content  of  57pe  in  magnetite,  and  the  reported 
weakness  of  the  resonance  signals,  it  appears  that  the  ZNR  method  has  very 
little  chance  of  success  as  a  means  of  measuring  rock  stress. 
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4.  ELECTRON  PARAMAGNETIC  RESONANCE  (EPR) 


This  type  of  resonance  (Bleaney  a  Stevens,  1953)  has  been  observed  in  ions  of 
three  groups  of  crystalline  salts  (among  other  substances),  as  follows: 

1.  Alums,  f luosi licates  and  Tutton  salts  of  the  "iron  group"  isotopes 
shown  in  table  7. 

2.  Rare  earth  salts  of  elements  from  atomic  number  57  (lanthanum)  to 
number  71  (lutecium).  These  are  in  actuality  not  particularly  rare  elements. 

3.  Actinides,  including  uranyl,  neptunyl  and  plutonyl  salts. 


Symbol 

Table  7. 

The  Iron  Group 

Isotopes  with  I  ^  0 

Natural 

abundance 

Atomic  number 

Name 

Mass  number 

Ti 

22 

Titanium 

47 

7.5% 

Ti 

22 

t! 

49 

5.5% 

V 

23 

Vanadium 

50 

0.24% 

V 

23 

It 

51 

99.8% 

Cr 

24 

Chromium 

53 

9.4% 

Mn 

25 

Manganese 

55 

100% 

Fe 

26 

Iron 

57 

2.2% 

Co 

27 

Cobalt 

59 

100% 

Ni 

28 

Nickel 

61 

1% 

Cu 

29 

Copper 

63 

69.% 

Cu 

29 

II 

65 

31.% 

EPR  is  possible  with  very  small  samples,  of  the  order  of  1  milligram,  because 
the  magnetic  moments  of  ions  are  about  1000  times  as  strong  as  those  of  nuclei. 
Hence,  there  appear?  to  be  a  possibility  that  the  larger  samples  (several  kilo¬ 
grams)  within  reach  of  the  magnetic  field  of  a  resonance  logging  tool  could  be 
effective  even  where  the  sensitive  isotope  is  only  a  trace  element  measured  in 
parts  per  million.  In  fact,  it  is  desirable  that  this  be  the  case,  since 
resonance  line  widths  must  be  as  narrow  as  possible  for  detection.  Murfunin 
(1965)  states  that  concentrations  must  lie  between  about  1  ppm  and  0.5  percent. 
Narrow  resonances  can  occur  only  when  each  paramagnetic  ion  is  surrounded  by 
several  layers  of  diamagnetic  atoms,  as  for  example  the  chromium  ion  Cr+3 
occurring  as  an  impurity  in  corundum,  AI2O3,  to  form  ruby.  The  iron  group 
elements  Ti,  Mn  and  Fe  are  among  the  more  abundant  constituents  of  accessory 
minerals  in  igneous  rocks.  Hence,  there  are  several  natural  minerals  capable 
of  EPR  response. 

Marfunm's  excellent  paper  (Marfunin,  1965)  summarized  all  known  work  on  EPR 
(ESR)  effects  in  natural  minerals  up  to  that  time,  including  Russian  measure¬ 
ments  on  over  200  mineral  species  from  numerous  localities.  An  enormous 
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literature  also  exists  on  ESR  of  numerous  atomic  species  dispersed  in  synthetic 
corundum,  fluorite,  periclase,  and  scheelite,  as  a  result  of  maser  and  laser 
studies  (Marfunin  and  Bershov,  1964).  Marfunin  found  that  at  room  temperature, 
the  great  majority  of  hi^  200  natural  specimens  showed  the  ESR  spectra  of 
Mn+2>  Fe+3  and  CR+3.  He  gives  tables  of  spin-llamiltonian  coefficients  deter¬ 
mined  by  ESR  for  several  common  rock-forming  minerals,  including  calcite, 
dolomite,  apatite,  coruncum,  rutile,  spinel,  fluorite,  and  beryl.  He  also 
observed  widespread  distribution  of  Fe+3  in  several  feldspars  and  in  certain 
rare  types  of  quartz.  These  data  would  require  further  study  for  practical 
application  to  our  problem. 

EPR,  like  NMR,  requires  a  steady  and  highly  uniform  magnetic  aligning  field, 
since  the  resonant  frequency  f  is  proportional  to  the  aligning  field  H,  as 
follows : 


f  =  gBH/2irh 


(7) 


Here  h  =  Planck's  constant  =  6.62  x  10~27  erg  sec 
B  =  Bohr  magneton  =  0.928  x  10“  erg/gauss 
g  =  spectroscopic  splitting  factor  =  1  to  2  +. 


EPR  has  customarily  been  performed  in  magnetic  fields  of  several  thousand 
gauss,  corresponding  to  microwave  frequencies,  for  several  practical  reasons. 

In  our  application,  however,  the  only  reasonably  uniform  magnetic  aligning 
field  available  is  that  of  the  earth,  which  is  of  the  order  of  only  0.5  gauss. 
This  corresponds  to  resonant  frequencies  of  0.7  Megahertz  (MHz)  to  perhaps 
1.5  MHz,  which  would  be  suitable  for  practical  stress  logging. 

The  EPR  line  width  is  determined  chiefly  by  the  ion  dilution  and  by  the  spin- 
lattice  relaxation  time,  which  is  a  function  of  the  mineral  and  the  temperature. 
For  this  reason  much  EPR  work  is  done  at  very  low  temperature,  which  is  not 
possible  in  our  application.  However,  through  a  fortunate  choice  of  mineral, 
it  is  conceivable  that  spin-lattice  relaxation  times  as  long  as  10-5  second 
at  room  temperature  could  be  found.  (Times  as  long  as  10_4  second  have  been 
achieved  at  very  low  temperatures.)  This  would  produce  a  line  width  of  0.1 
MHz,  which  is  compatible  with  EPR  frequencies  in  the  earth's  field.  At  about 
this  and  longer  times,  saturation  effects  begin  to  degrade  the  line  strength. 

The  mineral  chosen  must  also  contain  the  EPR  ion  in  sufficient  dilution  that 
the  spin-spin  relaxation  time  is  also  comparable  to  or  longer  than  10-5  sec¬ 
ond. 

Unfortunately,  EPR  in  many  substances  is  anisotropic;  the  g- factor  varies 
with  the  angles  between  the  crystal  axes,  the  r.f.  field  and  the  aligning 
field  H.  This  means  that  in  a  ploycrystalline  rock  where  all  possible  crystal 
orientations  occur,  and  where  the  direction  of  the  r.f.  field  varies  around 
the  sampling  coil  or  logging  tool,  the  resonance  will  be  smeared  out  and  its 
amplitude  reduced  for  magnetically  anisotropic,  paramagnetic  substances.  This 
further  restricts  the  choice  of  the  indicator  mineral. 

The  observability  of  EPR  effects  in  minerals  in  situ  will  also  require  ade¬ 
quate  resonance  amplitudes.  These  are  lower  at  low  frequencies  (Bleaney  5 
Stevens,  1953).  Hopefully,  this  effect  can  be  compensated  for  by  large  sample 
size. 
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We  did  not  plan  to  make  extensive  efforts  in  this  program  to  apply  the  EPR 

phenomenon  to  measuring  rock  stress.  The  outlook  was  not  considered  favorable 
for  the  following  reasons: 

1.  Tho  earth's  field  has  apparently  never  been  successfully  tried  to 
date  as  an  aligning  field  for  EPR. 

2.  The  choice  of  a  widely  useful  indicator  mineral  is  circumscribed 
by  constraints  on  mineral  composition  and  on  concentration  and  lattice  coup¬ 
ling  of  the  EPR-sensitive  species  that  are  still  largely  unknown. 


Nevertheless,  because  of  the  fact  that  trace  elements  have  produced  EPR  sig¬ 
nals  in  numerous  natural  minerals,  the  possibility  of  EPR  rock  stress  measure¬ 
ment  was  considered  open.  Near  the  conclusion  of  the  program,  as  it  began  to 
appear  that  strong  NQR  signals  could  not  be  obtained  from  common  minerals, 
further  study  of  EPR  was  performed. 


4.1  ELECTRON  FERROMAGNETIC  RESONANCE 

One  variety  of  EPR  which  initially  was  considered  of  interest  is  EFR  or  "ferro¬ 
magnetic  resonance"  (Toner,  1963).  Anti  ferromagnetic  or  ferromagnetic  reso¬ 
nance  has  been  observed  in  ilmenite  in  the  2  mm-4  mm  wavelength  range,  in  ruby, 
in  haematite  (microwave  and  infrared  regions) ,  in  synthetic  yttrium-iron  gar¬ 
net,  in  many  oxides  and  fluorides  of  manganese,  nickel,  chromium  and  cobalt, 
and  in  a  variety  of  alloys  (Praddaude,  1972).  In  general,  it  appears  that 
the  resonance  frequencies  are  too  high  for  our  application,  even  with  zero 
applied  magnetic  field  (figure  15).  The  reason  apparently  is  the  strong  in¬ 
teraction  between  the  unpaired-electron  spin  and  the  strong,  internal  magnetic 
fields  of  the  crystal  as  a  whole.  It  is  tentatively  considered  reasonable  not 
to  evaluate  further,  the  entire  class  of  minerals  exhibiting  EFR. 


4.2  EXPERIMENTS  ON  PARAMAGNETIC  MINERALS 

Nearly  the  entire  bulk  of  EPR  literature  titles  found  in  a  brief  search  of 
the  physics  and  chemistry  abstract  journals  were  concerned  with  synthetic  and 
organic  substances  not  occurring  in  nature.  Rather  than  perform  another  ex¬ 
haustive  literature  search  under  such  unpromising  conditions,  our  final  policy 
of  making  maximum  use  of  consultants  and  existing,  commercially-developed 
equipment  was  followed  once  again. 

A  Varian  model  4502-13  microwave  EPR  spectrometer  equipped  with  a  graphic 
chart  recorder,  a  large  magnet  and  stabilized,  adjustable-sweep  magnet  current 
supply  had  been  located  at  North  Texas  State  University  near  Dallas.  This 
instrument  operated  at  a  frequency  of  9000  MHz  and  was  equipped  with  a  tempera¬ 
ture-controlled  sample  cell  capable  of  operation  over  a  wide  temperature  range. 
Arrangements  were  made  for  Dr.  Paul  Jones  of  NTSU  to  test  numerous  samples  of 
selected  natural  minerals  containing  impurities  or  color,  in  this  equipment, 
lables  8  and  9  summarize  the  many  strong,  clear  resonances  recorded  in  the 
initial  tests.  Examples  of  the  recordings  are  given  in  figures  16  to  18. 

With  all  samples,  the  magnetic  field  was  swept  from  about  500  gauss  to  9500 
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Table  8.  EPR  "electron"  resonance  at  9000  MHz,  3300  gauss 


Estimated  voltage 
S/I  for  5  gm  sample 


Sample  No. 

Mineral 

and  Af  =  10  Hz 

Remarks 

21 

15 

Ottowa  sand 

Milky  quartz 

38  ( ?) 

30 

Others  also,  near 
2500  G 

5 

Smoky  quartz 

7 

I 

Orthoclase  feldspar 
(Fremont,  Calif.) 

60(?) 

Extrapolated 

9 

Orthoclase  feldspar 
(Goodsprings ,  Nev.) 

22 

3  other  large  reso¬ 
nances  at  1900-2700  1 

7 

Microcline  feldspar 
(Pringle,  S.  Dak.) 

13 

J 

Albite  feldspar 
(Gunnison,  Colo.) 

5.5 

10 

Hornblende 

8. 5(?) 

Another  near  2500  G? 

S 

Red  shale 

100  (?) 

Sample  conductive; 
0.05  gm  used 

Table  9.  EPR 

"manganese"  resonances  at  9000  MHz 

Sample  No. 

Mineral 

Estimated  voltage 

S/N  for  5  gm  sample 
and  Af  =  10  Hz 

Remarks 

R  Oil  shale 

(Colorado) 

14  Blue  calcite 

9  Brown  calcite 


500  7  strong  lines  be¬ 

tween  3100  5  3600  G 

188(?)  6  strong  lines  be¬ 

tween  2700  5  3000  G 
(misread?) 

40,000(?)  6  strong  lines  be¬ 

tween  3000  5  3600  G 


NOTE:  Symbol  (?)  indicates  a  large  extrapolation  based  on 
sensitivity  dial  settings. 
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Figure  16.  EPR  signals  from  brown  calcite,  gain  =  100,  x  -  0.005  sec 


gauss  over  a  5-minute  interval.  Any  resonances  seen  were  re-recorded  with  a 
narrower  magnetic  range  and  better  choices  of  gain,  filter  t  ard  sweep  time, 
to  obtain  good  resolution  of  the  resonance  width  or  slope. 

The  group  of  6  lines  centered  at  3,300  gauss  in  the  brown  calcite  (figure  16) 
can  probably  be  identified  by  their  number  and  field  value  as  due  to  manganese 
(Mn+2)  ion  impurity,  according  to  Dr.  Jones  (personal  communication);  the  six 
transition  energies  are  produced  by  splitting  of  the  electron  absorption  line 
by  coupling  to  its  nucleus  of  spin  I  =  5/2,  which  has  six  possible  states. 

This  identification  agrees  with  that  of  Weeks  (1973)  for  a  six-line  EPR  group 
found  in  a  feldspar.  The  single  resonance  found  near  3,300  gauss  (Lande  split¬ 
ting  factor  g  =  2)  in  most  of  the  minerals  can  be  attributed  to  single,  loosely- 
coupled  electrons  which  produce  color  in  many  compounds.  Figures  17  and  18 
show  this  resonance  well,  and  it  no  doubt  accounts  for  the  central  line  in  the 
group  of  7  in  oil  shale  (table  9). 

In  records  made  at  high  gain,  a  line  near  1525  gauss,  some  50  gauss  wide  always 
occurred.  This  was  found  also  with  an  empty  sample  tube,  and  is  attributed  to 
a  slight  contamination  by  iron  (Fe+^)  in  the  pyrex  glass  of  the  tube.  A  num¬ 
ber  of  other  signals,  including  broad  ones  such  as  that  seen  in  figure  18,  have 
not  yet  been  explained. 

Preparations  were  immediately  made  to  obtain  temperature  and  pressure  coef¬ 
ficients  of  some  of  the  observed  resonances.  Pressure  test  cells  were  made 
up  from  7  mm  o.d.  paper-phenolic,  epoxy-fiberglas  and  pyrex  glass  thick-walled 
tubing,  capable  of  being  placed  in  the  chuck  of  the  x-band  sample  cavity  of 
the  EPR  spectrometer.  Each  was  anchored  in  a  1/2-inch  pipe  bushing  by  threads 
or  a  flange  (on  the  glass)  and  sealed  with  epoxy  cement  or  a  neoprene  0-ring. 

A  compact  pressure  source  consisting  of  a  pair  of  high-pressure  needle  valves 
was  connected  to  the  cell  and  to  a  5000  psi  gauge,  with  the  entire  system  filled 
with  a  liquid  of  good  wetting  and  electrical  insulating  properties  (trichlor- 
ethylene) .  When  the  valve  stems  were  screwed  in,  their  displacement  readily 
developed  several  thousand  psi  of  hydrostatic  pressure  in  the  liquid.  The 
phenolic  tube  was  tested  to  4,700  psi  before  a  pinhole  leak  developed,  but  it 
proved  to  have  too  much  electrical  absorption  for  use  in  the  x-band  EPR  spec¬ 
trometer.  The  epoxy-fiberglas  tube  leaked  at  2,200  psi.  The  pyrex  tubes 
could  withstand  5000  psi  momentarily,  and  1000  psi  for  at  least  30  minutes, 
but  all  split  at  2000  to  3000  psi  after  about  5  minutes.  All  pressure  tests 
were  therefore  run  at  1000  psi,  with  occasional  extra  tests  at  2000  psi.  To 
prevent  contamination  of  the  microwave  cavity,  a  thin  gum-rubber  boot  was 
always  fitted  snugly  over  the  sample  tube.  This  boot  did  not  appear  to  pro¬ 
duce  any  extra  EPR  signals.  About  12  pyrex  sample  tubes  burst  within  this 
boot  without  mishap.  After  initial  tests  showed  that  the  valves  and  gauge 
were  subject  to  dangerous  magnetic  forces  when  attached  directly  to  the  sam¬ 
ple  tube  fitting,  a  2-ft  length  of  fine  stainless-steel  capillary  tubing  was 
used  for  interconnection.  Porous  plugs  of  paper  were  necessary  to  keep  all 
fittings  completely  filled  with  liquid  while  the  joints  were  made  up;  even  a 
small  air  bubole  would  prevent  the  pressure  from  developing.  Figure  19  shows 
the  final  pressure-testing  system  used. 

Table  10  summarizes  the  minerals  tested  in  the  pressure  and  temperature  cells. 

In  the  latter,  the  sample  was  alternately  cooled  to  30°C  and  electrically 
heated  to  200°C  by  a  stream  of  nitrogen  gas.  A  careful  comparison  of  the 
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Figure  19.  Arrangement  for  pressure-testing  of  FPR  frequencies 
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spectra  failed  to  show  any  temperature  shifts  of  frequency  greater  than  0.1% 
of  the  frequency  (or  field  strength) ,  which  was  the  limit  of  resolution.  How¬ 
ever,  weakening  and  broadening  of  some  lines  did  occur  at  the  higher  tempera¬ 
ture.  It  is  apparent  that  the  frequency  coefficients  of  temperature  for  all 
these  minerals  is  less  than  6  x  10_6/deg  C,  which  is  encouraging. 

Despite  careful  comparison  of  record  pairs  on  a  light  table,  no  frequency 
shifts  with  1000  psi  of  hydrostatic  pressure  could  be  seen  that  were  greater 
than  the  resolution  limit  of  0.3%,  except  for  a  -0.6%  shift  in  biotite  and  a 
+0.4%  shift  in  blue  salt  (which  had  received  10^  rads  of  gamma  radiation  in 
the  AEC's  Project  Salt  Vault).  If  verified,  these  shifts  are  very  substantial, 
of  the  order  of  3%  per  kilohar,  which  would  permit  usable  resolution  and  pre¬ 
cision  in  the  measurement  of  rock  stress.  We  were  unable  to  repeat  and  con¬ 
firm  these  measurements  for  lack  of  time  and  funds  under  this  program.  It 
would  be  desirable  to  re-run  all  the  FPR  pressure  tests  with  special  microwave 
fittings  to  which  much  greater  pressures  could  be  applied. 

Table  10.  Minerals  tested  for  EPR  shifts  with 
temperature  and  pressure 


Mineral 

Temperature 

Pressure 

Effects  seen 

Brown  calcite 

X 

X 

None 

Milky  quartz 

X 

X 

None 

Orthoclase 

X 

X 

None 

Hornblende 

X 

XX 

None 

Biotite 

- 

X 

-0.3%  shift 

Bituminous  coal 

- 

X 

None 

Blue  salt  (NaCl) 

( 1C^  rads  y-rays) 

- 

X 

+0.4%  shift 

Figure  20  shows  the  sharp  resonance  at  3296  gauss,  with  two  weak  satellites 
from  electron  interaction,  obtained  from  a  1-gm  sample  of  coal.  The  spacing 
of  these  satellites  should  be  pressure-sensitive,  but  they  could  not  be  seen 
during  the  pressure  experiment  because  only  a  1  mgm  sample  could  be  kept  sub¬ 
merged  in  the  trichlorethylene  in  the  1  mm  bore  of  the  pyrex  sample  tube. 
Figure  21  shows  the  complicated  spectrum  obtained  from  blue  salt;  those  with 
and  without  1000  psi  of  pressure  were  nearly  identical.  The  spectra  of  bio¬ 
tite  and  hornblende  were  qualitatively  similar  to  figure  IB. 
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4.3  DISCUSSION  ON  EPR  FOR  ROCK  STRESS  MEASUREMENT 


It  is  clear  that  strong  EPR  signals  can  be  obtained  from  many  minerals  with 
lab oratory -type  microwave  equipment.  There  are  indications  of  pressure  shifts 
of  the  resonances,  and  theoretical  grounds  for  expecting  pressure  effects,  at 
least  in  lines  arising  from  coupling  between  different  particles.  Investiga¬ 
tion  of  these  matters  in  the  literature  and  by  experiment  could  only  be  begun 
under  the  present  program. 

A  substantial  question  remains,  with  regard  to  the  electronic  measurement  of 
rock  stress,  whether  EPR  lines  can  be  obtained  and  the  effect  of  pressure  upon 
them  can  be  measured  at  frequencies  of  a  few  MHz,  using  only  the  earth's  mag¬ 
netic  field  for  alignment;  alternatively,  whether  a  synthetic  magnetic  field 
can  be  produced  around  a  borehole  by  properly  spaced  solenoids,  or  in  the  wall 
of  a  mine  by  a  large  loop,  that  will  be  homogeneous  enough  to  form  EPR  lines 
at  a  somewhat  higher  and  more  favorable  frequency.  One  or  the  other  will  be 
necessary  to  adapt  EPR  to  rock  stress  measurement. 

Some  guidance  has  been  obtained  from  the  literature  and  from  consultants. 
Bleaney  and  Stevens  (1953)  mention  that  many  EPR  line  widths  are  of  the  order 
of  6  gauss,  or  30  to  300  MHz,  so  that  a  somewhat  higher  frequency  is  required 
for  resolution.  Line  separations  are  also  substantial  (see  figure  16)  and  it 
may  be  necessary  to  present  two  lines  on  the  same  scan.  Furthermore,  signal 
amplitude  increases  with  frequency  (Orton,  1968),  partially  compensated  for 
by  the  accompanying  change  in  sample  size  (McMillan,  1968).  On  the  other 
hand,  microwaves  cannot  be  used  with  metallic  or  otherwise  lossy  samples,  and 
UHF  or  VHF  fields  must  be  used.  Decorps  and  Fric  (1972)  give  a  design  for  a 
sensitive  EPR  spectrometer  operating  at  40  to  210  MHz,  and  many  references  on 
the  subject.  Such  a  spectrometer  would  require  a  d.c.  magnetic  field  H0  of 
only  16  gauss  or  so! 

Inquiries,  particularly  to  Sam  Weiman  (of  the  University  of  Wisconsin  ?) 
brought  forth  a  number  of  encouraging  facts:  The  Varian  Company  once  produced 
a  60  MHz  spectrometer  employing  a  25  gauss  magnet  (see  Journal  of  Chemical 
Physics,  about  1963  -  author:  Acrivos) .  The  splitting  of  the  proton  line 
still  occurs  at  zero  field  in  some  materials,  and  gives  a  "sharp"  resonance 
at  54  MHz,  and  a  signal-to-noise  ratio  of  3  to  10  for  a  2  ml  sample.  Dr.  Dan 
Kivelson  at  UCLA  did  a  study  of  relaxation  processes  vs  pressure.  The  Philips 
Company  published  a  study  of  manganese  in  zinc  sulfide  at  UHF  frequencies 
before  1966.  James  Burgess'  thesis,  published  in  the  Physical  Re\iew  on  or¬ 
ganic  EPR  work,  employed  a  1000  MHz,  300  gauss  spectrometer  using  a  30  cm 
cavity  and  a  1-inch  by  10-inch  sample,  which  could  be  lossy.  We  have  not  been 
able  to  follow  up  these  promising  leads  under  the  present  program. 
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5.  RECOMMENDATIONS 


On  the  basis  of  this  preliminary  study,  the  Principal  Investigator  recommends 
further  R  and  D  on  the  possible  application  of  the  electron  paramagnetic  reso¬ 
nances  (EPR's)  found  in  most  impure  minerals,  to  the  measurement  of  rock  stress. 

Two  other  methods ,  NMR  of  hydrates  and  other  selected  minerals  in  the  earth's 
field,  and  NQR  with  very  large  samples  of  rock  containing  feldspars,  also 
appear  to  merit  some  additional  study. 

Further  investigations  of  NMR  in  ferromagnetic  accessory  minerals  such  as  mag¬ 
netite  (ZNR)  and  of  electron  ferromagnetic  resonances  (EFR)  are  not  recommended. 
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APPENDIX  1  to  TR  73-7 
PROPAGATION  OF  RF  FIELDS  TURO”3H  ROCK 


PROPAGATION  OF  RP  FIFLDS  THROUGH  ROCK 


Many  rocks  are  surprisingly  "transparent"  to  electromagnetic  radiation  in  the 
radio-frequency  range.  This  is  particularly  true  of  the  hard  rocks,  which  are 
of  special  concern  in  deep  mining  and  tunneling  projects  where  the  measurement 
of  rock  stress  is  important.  This  is  because  the  harder  rocks  (excepting  some 
porous  quartzites)  are  usually  massive  and  crystalline  and  contain  relatively 
little  uncombined  water.  The  moisture  content,  in  turn,  largely  governs  the 
electrical  conductivity  o,  which  in  turn  determines  the  rate  of  absorption  of 
a  propagating  electromagnetic  wave,  Kava  absorption  may  be  taken  as  a  measure 
of  the  propagation  loss  for  a  non-radi ative  (Coulomb)  r.f.  magnetic  field  from 
an  exploring  coil. 

The  attenuation  factor  a  determines  the  magnitude  vs  distance  x,  of  a  magnetic 
field  II  in  a  propagating  plane  wave  in  accordance  with  the  formulas: 

H  =  H0e-ax  ,  a  =(-)  y^J7  *  1207i(a/2)/T 

in  m  where  o  is  in  mhos/meter,  <  is  the  relative  dielectric  constant  and  u 
the  magnetic  permeability  is  assumed  to  be  close  to  that  of  free  space. 

There  is  an  r.f.  conductivity  o  equivalent  to  all  r.f.  losses  regardless  of 
mechanism,  which  can  readily  be  measured  on  an  r.f.  bridge.  A  fresh  rock  speci¬ 
men  is  cut  to  fit  as  the  dielectric  in  a  small  parallel-plate  condenser;  the 
capacitance  and  loss  tangent  are  measured  with  and  without  the  rock  dielectric, 
and  from  the  results  c(=i/p)  and  k  can  be  obtained  by  appropriate  formulas.  The 
process  is  fully  explained  elsewhere.*  From  these  results  a,  and  therefore  a 
typical  propagation  distance  D=6.9/a  for  60  dB  of  attenuation  can  be  calculated. 

Tables  A,  B,  C,  and  D  list  parameters  recently  measured  on  four  types  of  rock 
encountered  in  current  tunneling  projects: 

A  -  Massive  limestone  from  the  Gran  Sasso  vehicular  tunnel  project  near 
L'Aquila,  Italy. 

B  -  Granites  from  the  St.  Gothard  vehicular  tunnel  project  in  the  Swiss 

Alps. 

C  -  Schist  from  the  Washington,  D.C.,  subway  project. 

D  -  Red  shale  from  the  Irondequoit  sewer  tunnel  project  near  Rochester,  N.Y 

The  data  taken  at  1  MHz  are  of  greatest  interest  for  this  program.  Propagation 
distances  d  for  6  dB  of  field  attenuation,  a  reasonable  figure  for  stress  meas¬ 
urement,  can  be  obtained  from  the  right-hand  columns  of  the  three  tables  by 
means  of:  d  =  D/10.  For  example,  for  Sam' le  (f)  of  the  limestone,  soaked 
36  hours,  d  =  58.5  meters.  For  Sample  #66  z)  of  the  granite,  soaked  36  hours, 


Cock,  J.  C. ,  "RF  Flectrical  Properties  of  Bituminous  Coal  Samples," 
Geophysics ,  v.  35,  p.  1079-85  (1970). 

u (M> 
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Table  A  -  R.F.  Hlectrical  Properties  of  Gran  Sasso  Limestone 


Sample  (a) 
as  cut 

6 

2 ,  ohms 

K 

c  ,ohm  meters 

D/2  for  60  dB 

1  Ml!z 

24° 

10000 

7,1 

3050 

150  meters 

3  MHz 

15° 

2500 

5,0 

1220 

50 

25  MHz 

7,5° 

600 

4,0 

580 

22 

100  MHz 

1,5° 

164 

3,7 

790 

27 

Sample  (a) 

Soaked  4  hrs. 

100  MHz 

3,5° 

170 

3,6 

350 

12 

Sample  (a) 

Soaked  16  hrs. 

100  MHz 

5,5° 

165 

3,7 

220 

8 

Sample  (b) 
as  cut 

1  MHz 

13° 

15000 

5,0 

6000 

250 

5  MHz 

7° 

3500 

4,8 

2580 

100 

25  MHz 

4° 

750 

4,5 

960 

37 

100  MHz 

1° 

205 

4,1 

1060 

39 

Sample  (b) 

Soaked  4  hrs. 

100  MHz 

2° 

200 

4,2 

520 

19 

Sample  (b) 

Soaked  16  hrs 

100  MHz 

6° 

210 

3,4 

180 

6 

Sample  (c) 
as  cut 

1  MHz 

2,5° 

19300 

3,8 

46000 

1640 

5  MHz 

2,  ° 

3800 

3,5 

11400 

390 

25  MHz 

1,5° 

780 

3,6 

3100 

110 

100  MHz 

1°  (?) 

210 

3,4 

1260 

42 

Sample  (d) 
as  cut 

1  MHz 

3° 

21200 

4,1 

36500 

1350 

5  MHz 

2° 

4300 

3,7 

11000 

380 

25  MHz 

1° 

880 

3,15 

4500 

140 

100  MHz 

o 

Ln 

o 

•  o 

242 

3,35 

2500 

85 
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Sample  (e) 
as  cut 

6 

3  ,  ohms 

K 

,ohm  meters 

D/2  for  60  dli 

25  Nil  1  z 

2° 

880 

3 ,6 

2260 

80  meters 

ino  MHz 

1° 

200 

3,5 

1020 

35 

Sample  (f) 

as  cut 

1  MHz 

3,5° 

16300 

5 , 6 

23100 

1000 

5  MHz 

2,5° 

3230 

5,1 

6260 

259 

25  MHz 

2° 

690 

4,6 

1710 

67 

100  MHz 

2° 

170 

4,7 

424 

17 

Sample  (f) 

Soaked  4  hrs. 

1  MHz 

8° 

15000 

5,6 

8360 

363 

5  MHz 

4,5° 

3200 

5,1 

3570 

147 

25  MHz 

2° 

680 

4,7 

1700 

67 

100  MHz 

2° 

165 

4,8 

412 

16 

Sample  (f) 

Soaked  36  hrs. 

1  MHz 

5,5° 

15300 

5,4 

13800 

585 

5  MHz 

3° 

3180 

5,1 

5260 

217 

25  MHz 

n° 

L 

650 

4,2 

1620 

60 

100  MHz 

1,5° 

158 

5,2 

530 

7  ? 
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Table  B -  R.  F.  Electrical  Properties  of  St.  Gothard  Granites 


Sample  No.,  Frequency 

#66(a)  ,  wet  4  hrs  (IDO  MHz) 

Same,  soaked  15  hrs 

#67 (a)  ,  wet  4  hrs  (100  MHz) 

Same,  soaked  15  hrs 

#68 fa),  wet  4  hrs  (100  MHz) 

Same,  soaked  15  hrs 


6 

Z ,  ohms 

K 

p ,  ohm  meters 

D/2  f oi 

0,5° 

280 

2,4 

2200 

63 

O 

CM 

O 

V 

340 

2,6 

6700 

197 

0,5° 

295 

2,5 

2560 

74 

<0,2° 

290 

2,6 

6300 

186 

<0,2° 

290 

2,5 

6300 

182 

<0,2° 

300 

2,3 

6500 

180 

#66 (b)  ,  dry 


1  MHz 
5  MHz 
25  MHz 
100  MHz 


1° 

19000 

2,5 

144000 

4150 

1,2° 

3800 

2,7 

24000 

720 

0,5° 

750 

2,5 

14200 

410 

<0,2° 

185 

2,4 

6900 

195 

Same  ,  soaked  100  hrs 


1  MHz 
5  MHz 
25  MHz 
100  MHz 


15° 

16000 

3,0 

8200 

260 

7,3° 

3500 

2,7 

3500 

105 

2,5° 

725 

2,6 

2150 

63 

1,2° 

175 

2,7 

1100 

33 

#66 (c)  ,  dr) 


1  MHz 
5  MHz 
25  MHz 
100  MHz 


2° 

20500 

3,0 

61200 

1940 

1,5° 

4020 

3,08 

16200 

522 

0,7° 

820 

2,96 

7000 

220 

0,5° 

205 

2,96 

2450 

77 

Same,  soaked  36  hrs 

1  MHz 
5  MHz 
25  MHz 
100  MHz 


16° 

17200 

3,78 

6550 

5° 

3900 

3,20 

4700 

1,5° 

805 

3,06 

3220 

0,5° 

200 

3,10 

2400 

253 

154 

103 

77,5 
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Table  C  -  R.F.  Electrical  Properties  of  Washington,  D.C.  schist 


Sample  (a)  , 

tabular 

1.0  x  0.70  x 

0.20  inch,  f.f.  = 

0.65 

As  cut 

6 

Z,  ohms 

Y 

P ,ohm  meters 

D/2  for  60  dB 

1  MHz 

41.5° 

10,000 

7.4 

1,340 

67  meters 

5  MHz 

20.5° 

3,200 

5.4 

818 

35 

25  MHz 

8.0° 

760 

4.8 

485 

20 

100  Mlz 

4.0° 

190 

4.5 

245 

9.6 

Sample  (a). 

Soaked  45  hrs 

1  MHz 

37.5° 

11,000 

7.1 

1,610 

79 

5  mz 

18° 

3,250 

5.4 

940 

40 

25  MHz 

7.5° 

770 

4.8 

521 

21 

100  MHz 

4.0° 

192 

4.5 

245 

9.6 

Sample  (b) , 

tabular 

0.95  x  0.80  : 

c  0.20  inch,  f.f. 

=  0.71 

As  cut 

1  MHz 

19.5° 

14,000 

6.0 

4,050 

182 

5  Mlz 

10.5° 

3,500 

4.8 

1,850 

75 

25  NHz 

4.5° 

790 

4.2 

972 

37 

100  Mlz 

1.5° 

200 

3.9 

735 

27 

Sample  (b) , 
Soaked  45  hrs 


1  MHz 

34° 

11,100 

6.8 

1,920 

92 

5  MHz 

16.5° 

3,270 

5.  1 

1,110 

46 

25  NHz 

6.5° 

775 

4.4 

660 

26 

100  MHz 

3.0° 

200 

3.9 

368 

13 

Sample  (c) , 
As  cut 

wedge,  1 

.2  x  0.9  x  0.10-0.40 

inch,  f.f. 

=  0.80 

1  Mlz 

18.5° 

12,500 

4.8 

4,200 

167 

5  Mlz 

11.5° 

3,200 

3.9 

1,730 

62 

25  Miz 

5.5° 

730 

3.5 

820 

28 

100  NHz 

3.0° 

170 

3.6 

350 

12 

Sample  (c) 
Soaked  45  hrs 


1  Mlz 

44.5° 

7,400 

6.4 

1,140 

53 

5  Mlz 

24.° 

2,550 

4.9 

672 

27 

25  Mlz 

11.° 

640 

4.1 

361 

14 

100  MHz 

4.5° 

162 

3.8 

221 

7 

Notes : 

(1)  All  samples  alike;  grey  muscovite  schist,  foliation  minor  to  absent.  Very  tough. 
(21  Electric  field  perpendicular  to  cleavage  (foliation)  planes  in  all  samples. 

(3)  All  samples  dry-ground  to  size  on  SiC  wheel. 

(4)  Empty-cell  100  MHz  impedances:  a  5  b  =  437  0;  c  =  380  ft. 

(5)  Minimum  impedances,  obtained  with  wet  blotting-paper  pads  between  sample  and 
electrodes  and/or  with  5  psi  compression  on  foam  pads,  were  used  throughout. 
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Table  1)  -  R.F.  Electrical  Properties  of  Irondequoit  Red  Shale 


6 

Z,  ohms 

K 

p,  ohm-m 

D/2  for  60  dB 

Vert.  Polariz. 

25  MHz 

26° 

230 

8,3 

72 

3,8  meters 

100  MHz 

27,5° 

74 

5,5 

22 

0,95 

Horiz.  Polariz. 

25  MHz 

17° 

370 

8,15 

113 

5,9 

100  MHz 

21° 

100 

6 , 6 

25 

1,18 
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d  =  23.3  meters.  These  are  very  favorable  result  show  that  the  practical 

probing  depth  of  a  stress-measuring  r.f.  magnetic  field  in  such  rocks  would  be 
limited  only  by  field  spreading  and  coil  geometry. 

Measurements  were  not  made  at  1  MHz  in  the  fresh  red  shale,  but  the  25  MHz 
results,  d  =  0.38  m  and  d  =  0.59  m  can  probably  be  extrapolated  to  1  MHz  by 
multiplying  by  ten,  since  the  resistivities  are  comparable  to  those  of  grey 
shale  given  elsewhere.*  The  result  is  a  6  dB  propagation  distance  of  the 
order  of  3.5  meters  in  this  shale.  Even  this  would  be  satisfactory.  Shales 
contain  a  substantial  amount  of  finely  divided  feldspar,  which  might  make  them 
candidates  for  stress-measurements  through  the  nuclear  quadrupole  resonance  of 
aluminum-27. 


★ 

See  footnote  on  page  1. 
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